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CHAPTER  1 


INTRODUCTION 


With  increasing  space  activity,  the  use  of  light  and  flexible  structures  is  becoming 
important  to  reduce  the  high  cost  of  lifting  mass  to  orbit.  The  most  imponant  types  of 
structures  used  in  space  include  large  solar  arrays  and  antennas,  deployable  satellites  and 
structures,  precision  optical  systems,  and  space  stations.  Because  of  their  flexibility,  low  and 
closely-spaced  resonant  frequencies,  and  low  damping  coefficients,  advanced  structural  and 
control  techniques  are  required  to  achieve  the  stringent  pointing,  displacement  and  shape 
accuracy  requirements  for  space  structures.  It  has  been  found  that  passive  control,  combined 
with  active  control,  will  be  necessary  for  large  flexible  structures.  The  interrelationship  among 
the  various  disciplines  affecting  the  analysis  and  design  of  large  flexible  structures  is  presented 
in  Figure  1.1.  This  diagram  shows  that  the  aspects  of  modeling  of  the  structure,  control  system 
design  and  location  selection  of  control  devices  are  very  closely  related  and  play  central  roles 
in  the  flexible  structural  control  problem.  The  problem  of  control  system  design  of  large 
flexible  structures  can  be  broadly  divided  into  two  major  parts  -  modeling  and  controller 
design.  Flexible  structures  are  usually,  by  nature,  distributed  parameter  (continuous)  systems 
whose  dynamics  can  be  modeled  by  partial  differential  equations.  A  popular  method  to 
generate  a  finite-degree-of-freedom  model  for  approximating  a  flexible  structure,  which  has 
infinite  degrees  of  freedom,  is  the  finite  element  method  (FEM).  Because  of  the  large  number 
of  degrees  of  freedom  in  the  configuration  space,  the  finite  element  models  of  large  flexible 
structures  are  of  extremely  high  order,  typically  involving  thousands  of  variables.  It  is 
impractical  to  design  a  control  system  by  using  the  full  order  finite  element  model.  The  model 
reduction  procedure  is  used  to  reduce  the  high  order  finite  element  model  to  a  feasible  lower 
order  model  for  controller  design. 

The  model  reduction  procedures  might  introduce  the  observer  and  control  spillover 
problems  which  are  caused  by  the  unmodeled  modes  being  measured  and  excited  by  the 
controller.  The  instability  caused  in  a  simple  beam  because  of  spillover  was  demonstrated  by 
Balas  (1978).  Although  much  effort  has  been  expended  in  control  design  or  prcfilter  design  to 
filter  out  the  signal  from  the  unmodeled  modes,  the  fundamental  problem  has  not  been  totally 
solved.  Several  approaches  are  proposed  to  eliminate  the  effect  of  spillover  by  Longman 
(1979),  Sesak  and  Likens  (1979),  Lin  (1980),  Meirovitch  and  Baruh  (1981),  Czajkowski  and 
Preumont  (1987),  and  Preumont  (1988). 

As  stated,  the  problem  of  controller  design  is  a  major  issue  in  the  control  of  flexible 
structures.  We  can  classify  the  controller  design  into  two  major  types,  namely,  the  passive 
control  and  the  active  control.  The  passive  control  can  be  used  to  change  the  natural 
frequencies  and  damping  coefficients  only  by  small  amounts;  also  it  might  increase  the  stability 
margin  and  the  robustness  of  the  system.  The  active  control  can  modify  the  closed-loop 
performance  by  changing  the  damping  ratio  and/or  the  stiffness  to  some  extent.  However,  the 


active  control  performance  might  be  greatly  degraded  by  the  spillover  effect  of  the  residual 
modes  or  the  errors  and  uncertainties  present  in  the  system  model.  The  general  aspects  of 
structural  control  were  discussed  by  Leipholz  (1979)  and  Leipholz  and  Abdel-Rohman  (1985). 

To  meet  the  stringent  performance  requirements  of  flexible  space  structures,  several 
challenging  problems  should  be  overcome.  The  first  challenge  is  to  generate  an  appropriate 
reduced  order  model  for  representing  an  infinite  dimensional  distributed  structural  system. 
This  is  a  discrete  input  space-continuous  output  problem.  The  second  challenge  is  to  select  the 
number  and  locations  of  the  sensors/actuators.  This  selection  problem  directly  affects  the 
stability  of  the  closed-loop  system  (Balas  1978)  with  the  controller  designed  on  the  reduced 
model.  The  third  challenge  deals  with  robust  controller  design  which  retains  the  closed-loop 
stability  and  performance  requirements  under  some  amount  of  uncertainty  in  system 
parameters  or  design  variables. 

In  practice,  the  order  of  the  finite  element  model  involved  in  large  flexible  structures  is 
so  high  that  it  is  impossible  to  solve  a  control  design  problem  of  the  full  order  system.  Hence, 
the  model  reduction  problem  becomes  very  important.  To  describe  the  model  reduction 
problem,  consider  the  dynamics  problem  involved: 

x  =  Ax4-Bu  (1.1) 

y  =  Cx  (1.2) 

where  x,u,  and  y  are  the  state,  control,  and  the  output  vectors,  respectively.  A,  B,  and  C  are  the 
state,  input,  and  the  output  matrices,  and  the  dot  indicates  the  time  derivative.  From  this,  a 
reduced  order  model  (Ag,  Be,  Q)  with  reduced  order  state  vector  Xc  will  be  found  so  the  error 
criteria,  established  by  the  designer,  will  be  minimized. 

Xc=AeXc+BcU  (1.3) 

y  =  CcXc  (1.4) 

A  graphical  representation  of  the  reduced  order  model  is  shown  in  Figure  1.2. 

In  general,  the  model  reduction  methods  can  be  divided  into  two  main  categories  as: 
model  reduction  in  frequency  domain  and  model  reduction  in  time  domain.  Several 

comparative  studies  in  model  and  controller  reduction  were  made  by  Decoster  and 
Cauwenberghe  (1976),  Elrazaz  and  Sinha  (1981),  Hyland  (1984),  Sugimoto  et  al.  (1985), 
Hyland  and  Bernstein  (1985),  Parry  and  Venkayya  (1986),  Anderson  and  Liu  (1987),  and 
Oliver  (1987).  Most  of  the  existing  approaches  used  in  model  reduction  involve  an 
optimization  problem  which  finds  a  set  of  parameters  for  the  reduced  order  model  so  that  a 
stated  error  criterion  is  minimized.  The  model  decomposition  procedure  is  an  alternative  for 
the  modeling  of  large  flexible  structures.  The  substructures  technique,  which  has  been  used  in 
the  finite  element  method  for  solving  large  complex  structural  analysis  problems,  was  extended 


to  the  modeling  of  large  flexible  structures  for  control  design  by  Young  (1988)  and  Pan  et  al. 
(1989).  As  described  in  Pan  et  al.  (1989),  the  problems  of  model  reduction  and  the  spillover 
problems  are  avoided  by  using  the  substructures  technique. 

Since  10  years  ago,  the  literature  related  to  the  modeling  and  the  control  of  flexible  space 
structures  has  been  greatly  increasing.  Several  survey  papers  on  the  subject  have  been 
presented  in  the  past  by  Canavin  (1978),  Blair  (1978),  Dahlgren  and  Gunter  (1978),  Gran  and 
Rossi  (1979),  Seltzer  (1979,1980),  Balas  (1982),  Garibotti  (1984),  Nunre  et  al.  (1984),  and 
Santiago  et  al.  (1984).  The  aspects  of  passive  control,  decentralized  and  hierarchical  control, 
integrated  structural  and  control  design,  and  number  and/or  location  selection  of 
sensors/actuators  are  reviewed  in  the  following  sections. 

1.1  Passive  Control  in  Flexible  Structures 

The  passive  control  methods  for  vibration  suppression  have  been  studied  and  used 
extensively  for  a  long  time.  The  advantages  of  passive  control  are  as  follows: 

1.  It  is  simple  and  inherently  stable. 

2.  No  power-source  is  needed. 

3.  It  can  increase  the  stability  margin  of  the  overall  system. 

4.  The  passive  damping  augmentation  might  lead  to  reliable  and  robust  systems  with 
simplified  active  controls  requirements. 

5.  The  proper  design  of  the  passive  control  system  can  reduce  the  settling  time  of  the 
transient  response  and  the  peak-overshoot. 

Hughes  (1980)  and  Zak  (1988)  investigated  the  effect  of  passive  energy  dissipation  in 
large  space  structure  (LSS)  control  design.  In  the  passive  control,  addition  of  damping 
(damper)  and  stiffness  (spring  or  tendon)  is  usually  used.  Passive  damping  can  be  added  to 
flexible  structures  through  a  variety  of  approaches  including  the  addition  of  discrete  viscous 
dampers,  tuned-mass  dampers,  impact/fHction  joints,  and  constrained  layer  materials.  The 
applications  of  tuned-mass  dampers  in  high  rise  buildings  to  control  the  vibration  of  stmetures 
by  the  single-input  single-output  concept  were  discussed  by  McNamara  (1977),  Wiesner 
(1979),  and  Petersen  (1980).  Juang  (1984)  presented  a  technique  for  formulating  expressions 
of  the  optimal  tuning  law  for  an  elastic  system  including  a  truss  beam  and  a  tip  vibration 
absorber.  A  formula  relating  the  beam  parameters  to  the  size  of  the  tip  absorber  is  provided  for 
the  optimal  design  of  an  absorber  along  with  an  evaluation  of  its  performance  for  the  single 
mode  case.  Sesak  et  al.  (1986)  investigated  the  optimal  tuning  of  multiple  tuned-mass  dampers 
for  the  transient  vibration  damping  of  large  space  structures  by  using  modem  control  theory 
and  parameter  optimization  techniques.  Trudell  et  al.  (1980)  addressed  the  subject  of  passive 
damping  provided  by  the  use  of  viscoelastic  materials  for  large  space  structures.  In  the 
PACOSS  (Passive  and  Active  Control  of  Space  Structures)  program,  a  representative  system 
(Representative  System  Article)  was  developed  and  analyz^  by  Morgenthaler  and  Gehling 
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(1986).  The  modal  strain  energy  distribution  approach  was  used  to  determine  locations  where 
passive  damping  treatments  would  be  most  effective  in  a  structure.  Gehling  (1986)  also 
presented  the  benefits  of  passive  damping  with  regard  to  active  control  implementation  and 
retargeting  performance  of  a  large  space  structure. 

Damping  synthesis  for  flexible  space  structures  has  been  studied  by  Soni  and  Agrawal 
(1985),  and  Simonian  (1986).  Vegte  and  Hladun  (1973)  proposed  an  optimal  passive  beam 
vibration  control  design  by  using  optimal  control  techniques.  Prucz  (1987)  suggested  that  the 
passively  damped  joints  using  high  damping  viscoelastic  materials,  have  the  potential  of  being 
effective  practical  means  of  passive  vibration  control. 

A  review  of  the  damping  mechanisms  in  advanced  fiber-reinforced  composites  was  given 
by  Adams  (1986).  Prucz  and  Fu  (1988)  presented  the  feasibility  of  controlling  the  propagation 
of  mechanical  vibration  associated  with  elastic  stress  waves  in  fiber-reinforced  composites  by 
appropriate  tailoring  of  material  constituents  and  fiber  orientations.  Ashley  (1984)  presented  a 
study  of  the  size  effects  on  various  damping  mechanisms  and  thermal  dissipation  induced  by 
strain  gradients  during  vibration  of  monolithic  configurations.  Chen  and  Wada  (1988)  studied 
the  application  of  viscoelastic  damping  materials  employed  in  the  parallel  load  path. 

The  stiffness  modification  in  passive  control  is  usually  used  to  either  maintain  the  shape 
of  a  reflective  surface  (e.g.  Herbert  and  Bachtell,  1986)  or  reduce  the  vibration  amplitude  by 
increasing  the  stiffness.  Many  passive  tendon  control  applications  can  be  found  in  the 
structures  of  tall  buildings  and  long-span  bridges  to  achieve  higher  strengths.  Roorda  (1975) 
used  the  active  tendon  control  in  tall  structures.  An  experimental  tendon  control  system  for  a 
flexible  space  structure  was  investigated  by  Murotsu,  et  al.(1988).  Although  the  passive  control 
has  some  advantages,  it  is  usually  coupled  with  active  control,  to  obtain  what  is  known  as  a 
semi-active  or  hybrid  control,  to  achieve  the  stringent  dynamic  requirements  of  flexible 
structures. 

1.2  Decentralized  and  Hierarchical  Control 

The  models  of  large  scale  mechanical  and  structural  systems  are  usually  well  known  but 
often  very  complex.  The  system,  sometimes,  is  viewed  as  a  set  of  decoupled  subsystems 
(Siljak,  1979),  and  a  local  controller  (decentralized  control,  Singh,  1981)  is  synthesized  for 
each  subsystem  neglecting  the  coupling  among  the  subsystems  in  the  first  step.  To  satisfy  a 
specified  global  performance  requirement,  a  global  control  (hierarchical  control)  may  be 
introduced.  For  most  large  dynamic  systems  (such  as  LFS),  they  usually  can  be  described  as 
collections  of  N  interconnected  subsystems  (or  substructures).  For  the  i“  subsystem,  let  x,  be 
an  Ui  dimensional  state  vector,  Uj  be  an  mj  dimensional  control  vector  and  z^  be  an  q 
dimensional  vector  of  interactions  which  come  in  from  the  other  subsystems.  The  hierarchical 
structure  of  a  controlled  structure  is  presented  in  Figure  1.3.  The  subsystem  dynamics  are 
assumed  to  be  linear  and  can  be  represented  by  the  following  state  space  equations: 
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Xi(t)  =  Ai  Xi(t)  BiUi(t)  -1-  <^;Zi(t) 

(1.5) 

Xi(0)  =  Xio 

(1-6) 

It  is  assumed  that  the  interaction  vector  Zj(t)  is  a  linear  combination  of  the  states  of  the  N 
subsystems: 

j=i 

It  is  desired  to  choose  the  control  input  uj,  U2,...,  un  such  that  the  performance  index  J  is 
minimized.  The  performance  index  for  a  linear  time-invariant  system  can  be  defined  as: 


N 

J=I 

1/2  |xi(tf)|  -I-  l/2j[  Ixi(t)|  \  +  |ui(t)| 

^]dt 

i=l 

lo  ^ 

where  Xj  and  z;  satisfy  Eqs.(1.5)-(1.7),  and  |x|  =  x^  Qx,  and  Qi  and  Rj  are  positive 

semidefinite  and  positive  definite,  respectively.  Survey  papers  in  decentralized/hierarchical 
control  are  provided  by  Sandell  et  al.(1978),  and  Lindner  and  Riechard  (1986). 

A  heuristic  approach  for  the  control  of  serially  connected  dynamical  systems  with  and 
without  time  delays  between  the  subsystems  was  proposed  by  Singh  and  Coales  (1975).  The 
control  strategy  is  determined  by  minimizing  the  cost  function,  which  is  defined  as  the 
summation  of  output  error  squares  and  input  energy  for  each  subsystem  during  the  whole 
control  process.  Singh  (1975)  presented  a  two-level  algorithm  which  can  achieve  the  optimum 
decentralized  control  for  large  interconnected  systems.  A  comparison  of  two  hierarchical 
optimization  techniques  was  presented  by  Singh  and  Hassan  (1976).  Singh  et  al.  (1976) 
developed  a  decentralized  computational  procedure  to  find  the  optimal  feedback  gain  matrices 
of  high-order  linear  quadratic  problems.  The  interaction  prediction  principle  is  used  for  the 
development  of  the  closed-loop  decentralized  control.  Siljak  (1976)  proposed  a 
decomposition-aggregation  method  using  both  the  passive  and  the  active  stabilization  devices 
for  a  spinning  flexible  spacecraft  stabilization  problem.  A  decentralized  optimal  control  for 
stabilization  of  stabilizable  large-scale  systems  was  derived  by  Dceda  et  al.  (1983).  The 
optimal  control  in  a  decentralized  scheme  is  achieved  by  selecting  an  appropriate  performance 
index.  It  is  shown  that  optimality  of  the  control  is  preserved  for  a  modified  performance  index 
under  perturbations  in  interconnections  so  that  the  strength  of  coupling  does  not  increase. 
McClamroch  (1985)  presented  a  simple  form  of  hierarchical  control  in  the  vibration 
suppression  of  flexible  structures  including  the  dynamics  of  actuators.  First,  a  decentralized 
controller  for  the  flexible  structure  is  developed,  and  then  a  compensation  is  provided  for  the 
member  damper  actuators  to  suppress  the  effects  of  actuator  dynamics.  The  use  of 
decentralized  control  structures  for  LSS  was  investigated  by  Medanic  et  al.  (1987).  The 
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frequency  weighting  and  the  projective  control  techniques  are  used  to  ensure  that  the  control 
effect  does  not  spillover  into  the  unmodcled  dynamics.  Ozguncr  and  Yurkovich  (1987) 
proposed  an  approach  for  decentralized  frequency  shaping  for  incorporation  of  frequency- 
domain  (bandwidth)  constraints  into  the  control  design  of  a  large  flexible  space  structure  with 
the  dynamic  model  including  actuator  dynamics. 

The  decentralized/hierarchical  control  procedure  for  LFS  addressed  by  Janschek  and 
Surauer  (1987),  describes  a  multi-input  multi-output  design  procedure  based  on  a  well  known 
sequential  design  approach.  Local  Gow  level)  control  loops  for  baseline  stabilization  and  high 
level  control  loops  for  overall  performance  are  designed  to  satisfy  individual  criteria. 
Bernstein  (1987)  investigated  the  sequential  design  of  decentralized  controllers  by  using  an 
optimal  projection  approach. 

Ozguner  and  Yurkovich  (1986)  presented  a  control  strategy  involving  a  decentralized 
model  reference  adaptive  approach  (DMRAC)  using  a  variable  structure  control.  Experiments 
on  NASA’s  flexible  grid  were  performed  by  using  DMRAC.  The  DMRAC  method  was  also 
investigated  by  Lee  et  al.  (1988).  It  is  shown  that  this  method  can  achieve  cither  output 
regulation  or  output  tracking  with  adequate  convergence,  provided  the  reference  model  inputs 
and  their  time  derivatives  are  integrable,  bounded,  and  approach  zero  as  time  goes  to  infinity. 

Vukobratovic  and  Stokic  (1984)  presented  an  iterative  procedure  for  the  suboptimal 
synthesis  of  a  robust  decentralized  control  for  a  large-scale  nonlinear  mechanical  system  to  a 
set  of  nominal  trajectories  which  are  prescribed  together  with  a  region  of  allowable  parameter 
values.  A  decentralized  control  law  for  a  class  of  nonlinear  interconnected  systems  was 
proposed  by  Saberi  (1988).  This  optimal  decentralized  control  law  is  derived  by  defining  an 
appropriate  local  performance  index  for  each  isolated  subsystem  and  solving  these  local 
optimal  control  problems. 

West-Vukovich  et  al.  (1984)  presented  the  decentralized  robust  servomechanism 
problem  with  constant  disturbances/set  points  for  LFS.  It  is  shown  that  the  spillover  problem 
can  be  eliminated.  Lindner  (1985),  and  Davison  and  Gesing  (1987)  considered  the  control  of  a 
flexible  spacecraft  using  a  decentralized  technique.  The  model  order  reduction  in  decentralized 
control  was  discussed  by  Ozguner  and  Lee  (1983),  and  Yousuff  et  al.  (1986).  The  robustness 
property  of  reduced  order  decentralized  control  designs  was  studied  by  Young  and  Siljak 
(1985).  Additional  studies  on  the  decentralized  control  of  large  scale  systems  can  be  found  in 
Aoki  (1968),  Singh  and  Hassan  (1979),  Siljak  (1979),  Singh  (1981),  Xinogalas  et  al.  (1982), 
and  Davison  (1984). 

U  Integrated  Structural  and  Control  Design 

In  the  conventional  structural  control  area,  the  structural  design  and  the  control  design 
are  treated  as  two  separate  procedures.  Although  each  design  is  optimal,  based  on  the 
individual  criterion,  the  combined  system  might  not  be  optimal  in  the  global  sense.  Hence,  a 
great  deal  of  research  is  currently  in  progress  on  developing  methods  for  the  simultaneous  (or 
integrated)  design  of  the  structure  and  the  control  systems.  The  current  methodology  applied 
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to  the  integration  of  the  optimal  process  for  structures  and  controls  is  discussed  by  Wisshaar  et 
al.  (1986).  Cooper  et  al.  (1986A,1986B)  described  IMAT  (Integrated  Multidiciplinary 
Analysis  Tool)  for  the  integrated  control/structure  design. 

The  improvement  of  a  control  system’s  properties,  such  as  stability,  controllability  and 
sensitivity,  of  a  flexible  structure  by  changing  the  structural  parameters  has  been  studied  by 
Venkayya  et  al.  (1985),  Eastep  et  al.  (1986),  Haftka  et  al.  (1984,1987),  Khot  et  al.  (1984,1986), 
Flotow  (1986),  Schmit  (1987)  and  Rao  et  al.  (1989).  The  optimal  design  of  flexible  space 
structures  with  constraints  on  dynamic  properties  (or  response)  were  considered  by  Canfield  et 
al.  (1986),  and  Woo  (1987).  Some  optimization  techniques  used  in  structural  multiobjective 
design  were  discussed  by  Rao  (1984,1986,1987).  Yedav^li  and  Skelton  (1983),  and  Yedavalli 
(1984)  presented  the  determination  of  the  critical  parameters  of  LFS  with  uncertainty  by  using 
parameter  sensitivity  analysis  methods.  The  sensitivity  studies  in  controlled  structural  design 
were  presented  by  Gilbert  (1986),  Manning  et  al.  (1986),  and  Adamian  and  Gibson  (1987). 

The  weight  of  the  structure  was  minimized  with  constraints  on  the  distribution  of  the 
eigenvalues  and/or  damping  parameters  of  the  closed-loop  system  by  Khot  et  al.  (1985).  The 
structure/control  system  optimization  problem  was  formulated  by  Khot  et  al.  (1987)  with 
constraints  on  the  closed-loop  eigenvalue  distribution  and  the  minimum  Frobenious  norm  of 
the  control  gains.  A  unified  algorithm  for  sequential  (or  simultaneous)  design  modifications  of 
a  closed-loop  constant  gain  control  system  and  the  flexible  structure  to  be  controlled  was 
presented  by  Junkins  et  al.  (1984).  Hale  (1985)  considered  an  ellipsoidal  set-theoretical 
approach  to  the  integrated  structural/control  synthesis  for  vibration  regulation  of  flexible 
structures.  This  approach  attempts  to  maximize  the  allowable  magnitude  of  an  unknown  (but 
bounded)  disturbance  to  the  structure  while  explicitly  satisfying  specific  input  and  output 
constraints.  Both  structural  parameters  and  control  gains  are  variables  during  the  optimization 
process.  A  sequential  quadratic  programming  technique  was  used  to  solve  the  optimal 
structural/control  design  problem  by  Tseng  and  Arora  (1988).  They  solved  the  boundary  value 
problem  associated  with  the  optimally  controlled  structure  by  first  treating  it  as  an  initial  value 
problem  and  then  solving  it  by  using  the  iterative  method  of  nonlinear  quadratic  programming. 

Belvin  and  Park  (1988)  presented  a  method  for  the  optimization  of  the  closed-loop 
structural  system  using  only  structural  tailoring.  The  optimal  linear  quadratic  regulator  control 
theory,  in  conjunction  with  modal-space  control,  is  used  with  the  weighting  matrices  chosen, 
based  on  physical  considerations.  A  disturbance  rtKxiel  for  the  integrated  structural  and  control 
design  problem  was  proposed  by  Slater  (1988).  This  approach  uses  the  response  to  dynamic 
inputs  and  constraint  limits  to  establish  trade-offs  between  the  control  energy  and  the  structural 
mass.  Lim  and  Junkins  (1987)  presented  a  design  algorithm  with  numerical  applications  using 
stability  robusmess  measures.  Three  different  cost  functions,  namely,  the  total  mass  ,  stability 
robustness  and  the  eigenvalue  sensitivity,  have  been  optimized  with  respect  to  a  unified  set  of 
design  parameters  which  include  structural  and  control  parameters  and  actuator  locations. 

In  the  last  few  years,  the  simultaneous  control/strucnu’al  design  problem  has  been  treated 
as  a  multiobjective  optimization  problem.  It  can  be  formulated  as: 
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(1.9) 


subject  to 


Min  f(x,u) 

X.U 

g(x,  u)^0 


where  f  is  vector  of  objectives,  and  g,  x,  and  u  are  constraint,  structural  design,  and  control 
input  vectors,  respectively. 

Salama  et  al.  (198 

minimization  ,  in  structural  land  control  variables,  of  the  sum  of  the  structural  weight  and  the 
infinite  horizon  linear  regulator  quadratic  control  cost.  A  weighted-sum  method  for  solving  the 
multiobjective  optimization  problem  in  the  simultaneous  structural  and  control  design  of 
flexible  structures  was  proposed  by  Hale  et  al.  (1985).  A  weighted  total  mass  consisting  of 
structural  design  parameters  is  added  to  the  cost  functional  and  the  structural  parameters  are 
varied  to  find  the  minimum  total  cost,  and  the  mass  and  stiffness  distributions  of  the  structure 
are  determined  as  a  part  of  the  optimization  problem.  Onoda  and  Haftka  (1987)  presented  a 
similar  optimization  approach  to  solve  the  simultaneous  design  problem.  It  is  proposed  to 
minimize  the  weighted  structural  mass  and  weighted  input  energy  subject  to  the  output  norm 
constraints  by  finding  the  optimal  structural  parameters,  control  gains  and/or  weighting 
constants.  A  similar  problem  was  also  considered  by  Lust  and  Schmit  (1988).  Manning  and 
Schmit  (1987)  posed  the  integrated  design  problem  as  a  composite  objective  function  of  total 
weight  and  control  energy  subject  to  different  constraints.  The  composite  objective  function  is 
modified  to  a  single  objective  function,  such  as  control  energy  (or  total  weight),  with  a 
constraint  that  the  total  weight  (or  control  energy)  be  less  than  or  equal  to  the  specified  upper 
bound,  to  get  rid  of  the  scaling  problem. 

Rao  et  al.  (1988A,1988B)  investigated  the  design  of  actively  controlled  structures  using 
two  multiobjective  optimization  techniques:  goal  programming  and  game  theory.  The 
simultaneous  structural  and  control  design,  with  a  consideration  of  the  robustness  of  the 
controlled  structure  under  uncertainty  on  parameters,  was  proposed  by  Rao  et  al.  (1989).  In 
this  paper,  the  stability  and  performance  robustness  indices  are  introduced.  Several 
multiobjective  optimization  techniques  are  used  to  solve  the  three-objective  problem,  the  two 
robustness  indices  and  the  total  structural  weight,  to  get  a  more  robust  controlled  structure 
under  certain  physical  constraints.  However,  some  objective  functions  may  not  be  unimodal 
and  local  optima  will  be  obtained  by  conventional  optimization  techniques.  Hence 
revolutionary  optimization  techniques  need  to  be  developed  for  their  solution. 


6),  and  Miller  and  Shim  (1986)  considered  the  simultaneous 


S 


1.4  Sensor  and  Actuator  Location  Selection  Problem 


The  detennination  of  the  number  and  location  of  sensors  and  actuators  for  distributed 
flexible  structures  is  an  important  issue.  If  the  sensors  are  used  not  only  for  the  feedback 
control  gain  calculation,  but  also  for  monitoring  the  change  of  system  parameters  or  estimation, 
the  selection  of  the  number  and  location  of  the  sensor/actuator  becomes  more  critical. 
Generally  speaking,  this  problem  can  be  solved  using  optimization  techniques  or  cost 
decomposition  methods.  Using  optimization  methods,  the  problem  can  be  formulated  as: 

Min  f(n,m,ic,>^  (1.10) 

n,m,x.y 

subject  to  the  specified  constraints,  where  f  is  the  cost  function  which  is  usually  defined  as  the 
summation  of  the  output  error  covariance,  input  energy  and  the  cost  (weight,  budget)  of  sensors 
and  actuators.  The  variables  n  and  m  denote  the  number  of  sensors  and  actuators,  respectively. 
The  vectors  x  and  'y  represent  the  location  vectors  of  sensors  and  actuators,  respectively. 

A  survey  of  the  field  of  optimal  sensors  and/or  actuators  location  for  dynamical 
distributed  parameter  systems  modeled  by  partial  differential  equations  was  presented  by 
Kubrusly  and  Malebranche  (1985).  Omatu  et  al.  (1978)  solved  the  optimal  sensor  location 
problem  for  a  linear  distributed  parameter  system  by  using  a  criterion  to  minimize  the  trace  of 
the  optimal  filtering  error  covariance  function.  Schulz  and  Heimbold  (1983)  presented  an 
optimization  method  to  determine  the  locations  of  sensors,  actuators  and  feedback  gains  for  the 
control  of  flexible  structures.  Heuristic  integer  programming  was  used  to  select  the  actuator 
locations  in  large  space  structures  by  Haftka  et  al.  (1985).  Salama  et  al.  (1987)  cast  the 
location  selection  problem  as  a  combinatorial  optimization  problem  and  solved  it  by  an 
adaptation  of  the  simulated  annealing  heuristic  algorithm. 

Stieber  (1988)  investigated  the  interrelation  between  the  arrangement  of  sensors  and 
actuators  on  flexible  structures  and  the  hyperstability,  which  is  a  system  property  that 
guarantees  the  stability  of  interconnected  systems.  Velde  and  Caignan  (1984)  considered  the 
number  and  placement  of  control  system  components  including  possible  component  failures.  It 
is  solved  by  computing  the  performance  of  all  admissible  combinations  of  component 
locations.  The  placement  of  sensors  in  structural  control  and  the  associated  reliability  of  the 
system  was  analyzed  by  Baruh  and  Choe  (1988).  An  optimization  technique  is  used  to 
determine  the  best  locations  of  the  operational  sensors.  In  addition,  it  is  shown  that  modal 
filtering  by  means  of  spline  functions  is  more  desirable  for  modal  coordinate  extraction  than 
the  other  approaches  because  of  stability  and  robustness  considerations.  Montgomery  and 
Velde  (1985)  presented  the  sensor/actuator  location  problem  with  reliability  considerations. 
The  optimum  placement  of  sensors  and  actuators  for  static  deformations  and  shape  control  of 
LSS  was  investigated  by  Haftka  (1984),  and  Haftka  and  Adelman  (1987).  The  actuator/sensor 
failure  detection  methods  in  the  control  of  flexible  structures  were  studied  by  Baruh  (1986), 
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and  Baruh  and  Choe  (1987). 

The  component  cost  analysis  or  input/output  cost  analysis  was  developed  by  Chiu  and 
Skelton  (1981)  to  solve  the  optimal  number  and  location  selection  problem  in  LQG  systems. 
The  component  cost  at  each  possible  location  is  found  and  ranked  in  order.  The  possible 
locations  having  least  contribution  to  the  cost  function  are  truncated  and  the  remaining 
locations  are  considered  to  be  the  best  places  to  locate  the  sensors/actuators.  Additional  studies 
can  be  found  in  Chiu  and  Pan  (1983),  and  Skelton  and  DeLorenzo  (1983).  Lindberg  and 
Longman  (1984)  developed  a  method  to  determine  the  optimal  actuator  locations  for 
independent  modal  space  control.  A  projection  approach  was  developed  and  a  sensor  location 
criterion,  which  minimizes  the  error  caused  by  the  unobservable  subspace  under  nonstationary 
noise,  was  derived  by  Morari  and  O’Dowd  (1980).  Other  related  research  was  conducted  by 
Jai  (1986)  and  Lafontaine  and  Stieber  (1986).  Although  many  approaches  were  studied,  there 
is  no  guarantee  to  obtain  a  global  optimal  solution  for  the  sensor/actuator  location  selection 
problem  of  LFS. 

1.5  Outline  of  the  Report 

After  an  overview  of  the  recent  research  in  modeling  and  control  design  for  flexible 
structures,  the  report  is  organized  as  follows.  Chapter  2  introduces  the  interacting  substructure 
decentralized  control  design  approach.  The  idea  of  the  interacting  substructure  decentralized 
control  originated  from  the  substructure  techniques  used  in  structural  analysis  for  large 
complex  structural  systems.  The  interacting  force  balance  plays  the  key  role  in  the  individual 
control  gain  design  for  each  substructure.  The  control  gains  are  made  to  converge  after  several 
iterations  of  the  procedure.  This  method  does  not  need  model  reduction  and  hence  the  spillover 
problem  is  avoided.  A  passive  control  design  and  a  dual  passive/active  control  design  are 
introduced  in  Chapter  3.  In  Chapter  4,  the  stability  and  performance  robustness  indices  are 
defined,  and  the  integrated  structural/control  design  problem  is  considered  as  a  multiobjective 
optimization  problem.  The  optimal  actuator/sensor  location  selection  problem  is  addressed  in 
Chapter  5  and  three  methods  arc  considered  to  solve  the  problem.  Finally,  Chapter  6  offers 
some  concluding  remarks  on  this  research. 
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Figure  1.1  Relationship  Between  Techniques  of  Large  Flexible  Structures 
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Control  law  based 


reduced  order  model 


Reduced  Order  Model 


Residual  Dynamics 


Figure  1.2  Model  Reduction  and  Control  of  Large  Flexible  Structures 
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Figure  1.3 


Large  Rexible  Structure 


Hierarchicel/Decentreliied  Control  of  Large  Flexible  Structures 
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CHAPTER  2 


INTERACTING  SUBSTRUCTURE  DECENTRALIZED  CONTROL 


2.1  General  Formulation  of  Flexible  Structures 

For  flexible  structures,  the  equations  of  motion  can  be  expressed  as 

Mv  +  Cv  +  Kv  =  Du  (2.1) 

where  M,C,K,  and  D  are  the  mass,  damping,  stiffness,  and  the  input  matrices,  v  and  u  are  the 
displacement  and  the  input  vectors.  Letting  x*  =  {  v'  v  ),  Eq.(2. 1)  can  be  rewritten  as 

x  =  Ax  +  Bu  (2.2) 


where 


0  1 
-M“‘K  -m-'c 


(2.3) 


0 

M"’D 


(2.4) 


Without  loss  of  generality,  the  optimal  linear  quadratic  regulator  method  is  applied  to  design 
the  control  gain  of  the  feedback  controller  for  simplicity.  Thus,  the  input  vector  can  be 
presented  as 

u  =  -R"‘b'’'Px  (2.5) 

with  P  satisfying  the  following  matrix  Riccati  equation 

a’’'P+PA-PBR~’b'^P  +  Q  =  0  (2.6) 


where  Q  is  a  positive  semidefinite  output  weighting  matrix,  and  R  is  a  positive  definite  input 
weighting  matrix. 

2.2  Nodal  Condensation  for  Substructures 
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2.2.1  Direct  Condensation 


When  the  number  of  degrees  of  freedom  in  a  structure  becomes  very  large,  the  analysis  gets 
difficult  and  sometimes  becomes  impossible.  To  solve  this  kind  of  large  structures,  one  can 
divide  the  analytical  model  into  substructures  and  solve  the  smaller  problems  individually. 
The  equation  of  motion  for  a  specific  substructure  can  be  written  in  terms  of  the  internal  and 
the  boundary  degrees  of  freedom  as 


Maa  ^ab 

Xa 

+ 

f-aa  Cjb 

Xa 

+ 

Kaa 

Kab‘ 

Xa 

■fa‘ 

Mba  Mhb 

Xb 

Cba  Cbb 

Xb 

Kba 

Kbb 

Xb 

fb_ 

where  x  is  the  displacement  vector,  M,  C,  and  K  are  the  mass,  damping,  and  the  stiffness 
matrices,  and  f  is  the  force  vector.  The  vectors  and  the  matrices  in  Eq.  (2.7)  are  partitioned 
corresponding  to  the  boundary  degrees  of  freedom  (subscript  b)  and  the  internal  degrees  of 

freedom  (subscript  a).  Let  the  force  term  on  the  right  hand  side  of  equation  (2.7)  be  equal  to 
k  _  k 

so  that  the  system  response  can  be  expressed  as  ^^Vje^^*'.  For  each  force  component 
i=l  i=l 

there  is  a  corresponding  response,  thus  the  superposition  approach  can  be  applied.  The  original 
problem  then  can  be  divided  into  k  equations  of  motion  corresponding  to  any  specified  input 
frequency.  Considering,  for  simplicity,  only  one  force  component  f=fd“*,  Eq.(2.7)  can  be 
expressed  as 


Kaa+CaaO>-MaatO^ 

Kab+CabtO-MabCO^ 

Va 

'  " 

fa 

Kba+Cbati>-MbaW^ 

Kbb+CbbCO“Mbbt*>^ 

Vb 

fb 

Assume  that  no  force  is  applied  on  internal  nodes  (i.e.,  fg  =0)  so  that 

DaaVa  +  DabVb  =  0  (2.9) 

where 

Daa  =  Kaa  +  C^(0  "  MaaCO^  (2. 10) 

Dab  =  Kab  +  CabtO  -  MabW^  (2.11) 


The  exact  transformation  from  the  degrees  of  freedom  at  internal  nodes  to  those  of  boundary 
nodes  can  be  obtained  by  solving  Eq.(2.9). 
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Vj  —  D  aal^ab^b 

(2.12) 

Let  O  be  the  fixed  interface  modal  matrix  and  be  a  diagonal  matrix  consisting  of  the  natural 
frequencies  corresponding  to  the  modal  matrix.  The  following  relations  between  modal  data 
and  mass,  damping,  and  stiffness  matrices  are  satisfied. 

<D^Maa<l)  =  Identity  Matrix 

(2.13) 

(2.14) 

<l)'^C,a<I>  =  2'P 

(2.15) 

where  'P  is  damping  coefficient  matrix.  The  damping  is  assumed  to  be  proportional  damping  so 
that  'P  is  diagonal.  Substituting  these  relations  into  Eq.(2.10)  yields 

D-‘  aa  =  0(i2^  +  24^0)  - (I>T  (2. 1 6) 


In  Eq.(2.16),  the  term  in  parenthesis  is  in  diagonal  form  by  assuming  that  all  eigenvalues  are 
distinct  and  hence  the  inverse  of  it  is  the  inverse  of  each  diagonal  term.  For  example,  the  i*** 
term  can  be  expressed  as 

(Q^  +  2'Pia)-a)^r*  = 


2'Pi 

1 - + 

4H>f  1 

2 

0)^  - 

£2^ 

>  < 

- 

(Q2  +  2'P(o  - 

00 

1 

4'P? 

0) 

3 

0)^  ► 

£2f 

£2f 

£2f 

(2.17) 


Substituting  Eq.(2.17)  into  Eq.(2.16),  the  inverse  of  transfer  matrix  Dg*  can  be  obtained  as 
follows. 


(2.18) 


Dli=K;i 


1-Ca.K 


■^(0+  [MaaKli  +(CaaKli)2j 


0)‘ 


+  <I>  diag  [residual  part]<l>  co 


where  the  residual  part  is  contained  in  the  last  term  of  Eq.(2.17).  Thus  the  transformation 
matrix  can  be  written  as 


Va 

-DL’Dab' 

Vb 

^  J 

I 

L  J 

(2.19) 


This  transformation  is  a  function  of  the  frequency  of  excitation.  The  frequency  content  of 
excitation  is  usually  unknown  in  advance.  If  the  frequency  content  of  excitation  is  small 
compared  to  the  system  natural  frequencies,  the  frequency  dependent  terms  in  the  transfer 
matrix  of  Eq.(2.19)  will  be  less  important  and  a  simplified  transformation  can  be  obtained  as 
follows. 

D-'  =  (2.20) 


Dab  ~  ^ab 


Va' 

-K^ia'Kab' 

Vb 

N  J 

I 

k.  ✓ 

Vb 


(2.21) 


(2.22) 


As  a  more  general  case,  a  serial  structure  is  considered  with  the  equation  of  motion  of  each 
substructure  in  a  tridiagonal  form  as 


VBk 

CBkBk  CB^Ak  0 

VBk 

^AtBk  ^AiAk  ^AkBk„ 

^Ak 

+ 

^AkBk  ^A^Ak  ^AkBk,i 

M 

< 

■> 

0  MB^„Ak  MB,^.Bk.. 

VBk,. 

0  ^Bk,,Ak  ^Bk,,Bk,i 

^Bk,. 
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KBkBk 

KBkAk 

0 

VBk 

Fflk 

+ 

^AkBk 

KAkAk 

*^AkBk., 

< 

> 

= 

FAk 

0 

Kfik^Ak 

^Bk,iBk*i 

VBk,. 

FBk., 

in  which  the  subscript  denotes  the  internal  degrees  of  freedom  of  substructure  k,  and  Bk  and 
Bk+i  denote  the  boundary  degrees  of  freedom  between  substructure  k  and  connected 
substructures  k-I  and  k+1,  respectively.  For  a  single  frequency  excitation,  Eq.(2.23)  can  be 
rewritten  as 


DfikBk 

DfikAk  0 

''Bk  ' 

‘‘Bk 

l^AkBk 

DAkAk  DAkBk., 

VAk 

= 

fAk 

0 

Dfik.iAk  DBk.iBk.i 

>•> 

fBk.. 

(2.24) 


The  internal  degrees  of  freedom  can  be  reduced  to  boundary  degrees  of  freedom  by  assuming 
equal  to  zero: 

VAk  =  -D"'  AkAk  (DAkBk  Vfik  +  DAkBk,,  ^Bk,,  )  (2-25) 


Thus  the  transformation  can  be  obtained  as 


''Bk 

''Ak 

''Bk., 

s 


iBk  0 

-D"‘AkAkDAkBk  -D~*AkAkDAkBkM 

0  iBk,. 


(2.26) 


where 

I^AkAi  =  ^A^Ak  +  CAkAk^B-  Ma^Ai^^ 


DAkBk  =  f^AkBk  ■••CAkBk®"  ^^AiBk®^ 


(2.27) 


^AkBk,.  =  KAkBk.,  +  CAkBk.,  W  -  MAkBk,, 

The  transformation  matrix  can  be  reduced  to  a  simple  form  if  the  frequency  content  of 
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excitation  is  small  compared  to  the  system  natural  frequencies.  Thus 


Tk  = 


Ib,  0 

-K“'AkAkKAkBk 

0  Ib.. 


(2.28) 


2  J.2  Indirect  Condensation 


An  alternative  condensation  approach,  using  modal  data  of  each  substructure,  is  developed 
in  this  section.  The  full  order  modal  matrix  (<1>)  need  to  be  found  in  this  approach  and  the 
modal  coordinates  arc  denoted  by  the  vector  y.  The  displacement  v  can  then  be  represented  as 

v  =  <I>y  (2.29) 


and  the  inverse  of  exists,  so  that 


y  =  <I>~*v 


(2.30) 


Suppose  that  the  modal  coordinate  vector  y  is  arranged  in  two  parts  which  are  truncated  modes 
(y,)  and  retaining  modes  (y,).  The  physical  coordinates  (v)  are  also  arranged  in  two  parts 
which  arc  internal  degrces  of  freedom  (Va)  and  boundary  degrees  of  freedom  (Vb).  Then 
Eq.(2.30)  can  be  rewritten  in  the  following  form 


<*>ta 

^rb 


(2.31) 


where 


a)-‘  = 


<t>ta  <I>lb 
^ra  ^rb 


(2.32) 


In  Eq.(2.3I),  the  number  of  truncated  modes  will  be  chosen  same  as  the  number  of  internal 
degrees  of  freedom.  To  eliminate  the  truncated  modes,  vector  yt  is  set  to  be  zero  and  hence 


<*>laVa  +  <I>ibVb  =  0 


(2.33) 


The  relation  between  Vg  and  Vb  then  can  be  obtained  as 

Va=-^  *ia<I>ibVb  (2.34) 

and 
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(2.35) 


where 

Tab  =  u  (2.36) 

23  Interacting  Substructure  Decentralized  Control  for  Serial  LFS 

For  a  large-scale  interconnected  system  composed  of  a  number  of  subsystems  or  for  a 
system  distributed  widely  in  space,  the  scheme  of  decentralized  control  has  advantages  in 
computation  and  implementation  of  control  laws  (Siljak,  1976,  Sandelljr.  et  al.,  1978,  and 
Singh,  1981).  When  decentralized  optimal  control  is  used,  the  resulting  closed-loop  system  is 
guaranteed  to  have  robust  stability  properties  against  variations  in  the  open-loop  dynamics 
(Ikeda  et  al.,  1983).  Hence  the  stabilization  of  large-scale  systems  by  means  of  local  state 
feedback  has  been  applied  extensively. 

In  this  section,  a  new  decentralized  control  approach  for  serial  LFS  is  developed  by  using 
interacting  substructure  technique.  The  optimal  linear  quadratic  regulator  method  is  used  for 
the  controller  design.  Let  the  complete  flexible  structure  be  divided  into  n^  substructures  as 
indicated  in  Fig.  2. 1  and  the  equations  of  motion  of  the  k^**  substructure  be  given  by  Eq.(2.7) 
which  can  be  rewritten  in  state  form  as  in  Eq.(2.2).  The  conttoller  design  starts  from  the  first 
substructure  which  contains  nodes  belonging  to  the  global  boundary  as  well  as  the  local 
boundary  between  substructures.  Without  loss  of  generality,  optimal  LQR  method  is  used  in 
control  design  for  simplicity.  In  the  first  iteration,  the  controller  is  designed  by  assuming  the 
local  boundary  to  be  free.  The  siate  equation  and  the  control  law  can  be  written  as  follows: 

xi  =  AjXi -t- BjUi  (2.37) 

The  optimal  state  feedback  gain  is  given  by 

ui  =  -G,  xi  =  -RT‘  b|  P,  xi  (2.38) 

where  Pi  satisfies  the  matrix  Riccati  equation 

A|Pi  -t-  P,  A,  -  Pi  B,  R?'  B|Pi  +  Qi  =  0  (2.39) 

The  state  vector  xi  can  be  divided  into  interior  (xi^)  and  boundary  (xig)  states  so  that, 

,.T  .T  .,T  ..T  T 

Xi  =(xa,  XB,  xa,  XB,)  =AiXi-i-BiUi 
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=  Ajxi  + 


0 

-BiG, 


Xl 


(2.40) 


where  Bi  is  lower  half  part  of  matrix  Bj.  The  control  force  can  be  obtained  by  multiplying 
-BiGi  by  Ml  and  expressing  in  the  following  form: 


fA. 


=  -MiBiGiXi 


j 

(2.41)  ; 


! 


where  fA,  and  fe,  are  equivalent  forces  applied  on  internal  and  boundary  degrees  of  freedom. 
Applying  the  transformation  in  Eq.(2.22)  gives 

rTA,B.  0 


r  -N 


‘A. 

^B, 


=  — Ml  BjGi 


I  0 

0  Ta,b. 
0  I 


XB. 

XB, 


(2.42) 


where 


Ta,b,=-K-'a,a.Ka,b,  (2.43) 

The  equivalent  forces  applied  on  boundary  degrees  of  freedom  can  be  expressed  as  Fb 


fA. 

Fb  = 

OIb,  j 

ffl, 

»  * 

=  flB,XB,  +f2B,XB, 


(2.44) 


where  Ib,  is  the  identity  matrix  with  the  same  dimension  of  Xb,  . 

The  intermediate  substructure  k  (k^l  and  nj)  will  have  boundaries  connected  to  the 

(k-l)‘**and  the  (k+l)***  substructures.  Let  the  corresponding  boundary  states  be  represented  as 

xb^  and  xb^.,  and  the  boundary  force  contributions  be  denoted  as  Fb^  and  Fb^,,  ,  respectively. 

The  equations  of  motion  of  the  k*  substructure  are  given  by 

r.T  .T  .T  ..T  ..T  ..T  IT 
Xk  =  [XBk  Xa^  Xb^.,  Xb^  Xa^  Xb^.,  J 
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0 

0 

=  AfcXk  +  BkUk  +  1 

Mt'FB, 

+ 

(2.45) 


In  the  first  iteration,  Fb^,,  (k=l  to  Os-l)  is  assumed  to  be  a  zero  matrix  so  that  the  boundary 
forces  can  be  explicitly  expressed  as 


and 


+f2BkXBj 

0 

0 


(2.46) 


0 

0 

flB».,XB,.,  +f2Bk„XB, 


(2.47) 


Substituting  Eqs.(2.46)  and  (2.47)  into  Eq.(2.45)  gives 

Xk  =  (Ak  +  Fb^  +  Fb^^,  )  Xk  +  BfcUk 


=  AkXk  +  BkUk 


where 


FB,Xk  = 


0 

Mk’FB, 


(2.48) 


(2.49) 


FB,.,Xk  = 


0 


(2.50) 


The  optimal  linear  quadratic  regulator  method  is  applied  to  solve  the  constant  optimal  feedback 
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gain  for  substructure  k: 


ujc  -  -GjfXk  -  -Rj('  B J P^xic 


(2.51) 


where  satisfies  the  Riccuti  equation 

AfcPk  +  PkAic  -  PkB|cRk*BkPic  +  Ok  =0  (2.52) 


Similarly,  the  forces  applied  on  the  boundary  which  are  contributed  by  the  actuators  involved 
in  the  substructure  k  can  be  formed  as 


BfcUk  = 


r 

0 

-BkGk 


(2.53) 


where  Bk  is  the  lower  half  part  of  matrix  Bk-  Multiplying  -BkGk  by  the  mass  matrix  gives  the 
control  force  vector; 


Fb. 

Fa, 

Fb,., 


--MkBkGkXk  =-MkBkGkTk 


XB, 

XB, 


(2.54) 


where  Tk  can  be  obtained  from  Eq.(2,28)  for  direct  condensation  or  Eqs.(2.35)  and  (2.36)  for 
indirect  condensation  as: 


Tk  0 
0  Tk 


(2.55) 


The  forces  acting  on  the  boundary  degrees  of  freedom  are  given  by  Fb^  and  Fb^^,  in  Eq.(2.54). 
In  order  to  include  boundary  forces  into  the  equations  of  motion  of  connected  substructures,  the 
equivalent  boundary  forces  should  be  transformed  to  be  functions  which  depend  on  its  own 
boundary  degrees  of  freedom.  From  direct  condensation  Eq.(2.24),  the  first  and  the  third 
equations  give 
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^BkHk  +  ^BkAk  VAi  =  Ffik  (2.56) 

DBk.,  Ak  VAk  +  DBk.,Bk,,  VBk,,  =  FBk„  (2.57) 

Substituting  VAk  from  Eq.(2.25)  gives 

(DBkBk  +  De^AkTAkBk)  VBk  +  DfikAkTAkBkM  ^Bk.,  =  ^Bk  (2.58) 

OBk.,  AkTAkBk  VBk  +  (DBk.iBk.,  +  DB,.,AkTAkBk.,  )  ^Bk,,  =  Ffik.,  (2  59) 


where 

TAkBk=-D-‘AkAkDAkBk  (2.60) 


Assuming  Fb^  or  Fsk^,  to  be  zero  gives  the  relationship  between  the  displacement  vectors  vr^ 

andvBk,,: 

VBk  =-(DBkBk  +E>BkAk'FAkBkr‘  DBkAkTAkBk*, ''Bk.,  = ’^BkBk.,  VBk„  (2.61) 


'^Bk.,  =-(DBk.,Bk.,  +F>Bk.,AkTAkBk.,)  '  DBk^AkTAkBk ^Bk  =TBk.,BkVBk 


(2.62) 


The  transfer  matrices  in  Eqs.(2.61)  and  (2.62)  can  be  reduced  to  the  corresponding  stiffness 
matrices  if  excitation  contains  only  low  (compared  to  the  system  natural  frequency)  frequency 
components  (i.e.,  D  is  substituted  by  K  in  Eqs.(2.61)  and  (2.62)).  Thus  the  equivalent 
boundary  forces  can  be  rewritten  as 


pBk  =  ~  ^  MitB|cG|cTkT|(+|j( 


(2.63) 


Fb 


k*l 


»■  ' 

’‘Bk., 

0  01b,.,  MkBkGkTkTk.k^, 

L  J 

’‘Bk., 

k  * 

(2.64) 
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where 


Ib,  0 

Tb,..b,  0 

0  Ib, 

0 


(2.65) 


and 


Tk,k+1 


Tb^b^.,  0 

^B,.,  0 

0  Tb^b^., 

« 


(2.66) 


Equations  (2.63)  and  (2.64)  give  two  equivalent  forces  which  are  contributed  by  the  actuators 
of  substructure  k  and  act  on  the  boundary  nodes  of  substructures  (k-1)  and  (k+1),  respectively. 
The  transformation  matrix  for  indirect  condensation  approach  is  much  simpler.  For  the 
interacting  force  on  boundary  Bk+i,  only  the  degrees  of  freedom  belonging  to  are 
considered  as  boundary  d.o.f.’s.  Hence,  the  transformation  matrix  in  Eqs.(2.35)  and  (2.36)  are 
valid  after  proper  assignment  of  interior  and  boundary  degrees  of  freedom.  The  controller  of 
each  substructure  can  be  designed  one  after  the  other  sequentially.  In  the  first  iteration,  the 
acting  force  contributed  by  the  connected  substructures  is  on  one  side  only  (except  the  first 
substructure).  In  the  consecutive  iterations,  all  substructures,  except  the  first  and  the  last,  are 
subject  to  interacting  forces  on  both  sides.  The  control  gain  of  each  substructure  converges 
iteratively.  Proper  summation  then  can  be  used  to  obtain  the  global  control  gain  corresponding 
to  the  complete  system  state  Xj  which  is  the  union  of  Xi,X2,...,Xn^. 

2.4  ISDC  for  Substructures  Containing  No  Actuators 

The  basic  idea  of  ISDC  discussed  in  the  last  section  is  the  interacting  force  balancing.  If 
there  is  any  substructure  containing  no  actuator,  the  interacting  force  to  the  connected 
substructures  is  assumed  to  be  zero  which  can  not  balance  with  those  forces  from  connected 
substructures.  Hence  the  procedure  described  in  the  last  section  is  not  valid  any  more.  To 
generalize  the  ISDC  procedure,  interface  stiffness  and  inenia  loading,  used  in  the  dynamic 
analysis  of  large  complex  structures,  is  used  to  connect  the  substructures  containing  no 
actuators.  For  example,  two  connected  substructures  1  and  2  have  the  same  boundaiy  (i.e., 
vb,=vbj)  and  the  equation  of  motion  (EOM)  of  substructure  2  is  influenced  by  substructure  1. 
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Thus,  the  modified  EOM  of  substructure  2  can  be  written  as: 


MajAj 

MajRj 

CajAj  ^AjB, 

> 

1 

+ 

MbjAj 

MbjBj+Mb, 

VB, 

CbjAi  CbjBj+Cr, 

KajAi  KajB, 

va, 

fA.' 

KbjAj  KrjBj+Kb, 

VB, 

fB, 

(2.67) 


where 


Mb,  =T'fi2M,T,2 

(2.68) 

Cb,  =T^  12C1T12 

(2.69) 

Kb,  =tT,2K,T,2 


(2.70) 


and  the  transformation  matrix  T12  can  be  obtained  from  Eq,(2.22)  as  follows: 

r-(KA,A,r’KA.B. 


T,2  = 


I 


(2.71) 


where  Mj ,  Ci ,  and  Ki  are  the  mass,  damping  and  the  stiffness  matrices  of  substructure  1 . 
For  a  serial  structure,  a  more  general  form  can  be  derived  for  substructure  k: 


— k-l 

' 

MBkBk+M  Bk  MflkAk  0 

VBk 

MAkBk  MA,Ak  MAkBk,, 

VAk 

— k+1 

0  MB^  ,Ak  '^Bk.,Bk.,+M  Bk., 

VBk,. 
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CbjBi+C  Bk  CB^Ak  ^ 

+  ^AkBk  ^AkAk  ^AkBk.i 

— k+1 

^  ^Bk.kAk  CBk..Bk.,+C  Bk., 

KBkBk+K  Ke^Ak  0  vb^  Fb^ 

+  KA,Bk  KA,Ak  KAkBk.,  VAk  =  pAk  (2.72) 

^Bk,,Ak  KBk„Bk*,+K  Bk„ 

where 


M  Bk  =T^^k-l.kMk-lTk_i,k 

(2.73) 

^  Bk.i  ^'T^k+l.k+lMk+lTk+i.k+i 

(2.74) 

C  Bk  =  k-l.kCk-l^k-l.k 

(2.75) 

C  Bk.i  =T^k+l.k+lCk+lT’k+l.k+l 

(2.76) 

K  Bk  =T^  k-l.kKk-lTk_i,k 

QJl) 

— k+1  j 

K  Bk.i  =T  k+i.k+lKk+iTk+i.k+i 

(2.78) 

Transformation  matrices  Tk_i.k  and  Tk+i,it+i  can  be  derived  from  Eqs.(2.28),  (2.61),  and  (2.62) 
and  expressed  as 


27 


T'k-l.k  - 


I 


-(KA,,A,,r'KA,,B 


‘A,.,AnKA^.,nk.,)  '  KBk.,Aj.,  K  'Ai^.A^.,  Ka^  ,0^ 


r-\ 


(2.79) 


Tk+i,k+i  - 


(^Bk,2Bk,i  ^Bk.jAk*!  ^  Ak,iAk»i  ^Ak.iBk.i^  ^Bk,2Ak,i  ^  Ak.iAk.i  ^Ak.iBk.j 

AktiAk.i  ^^Ak,iBk,i 


(2.80) 


Thus  the  substructures  which  contain  no  actuators  can  be  condensed  to  the  connected 
substructures  using  the  equivalent  stiffness  and  inenial  loading.  The  interacting  force  used  in 
ISDC  approach  can  then  be  transmitted  through  the  substructure  containing  no  actuator  to  the 
other  substructures.  For  example,  let  the  k***  substructure  be  the  only  substructure  containing  no 
actuator.  Then  the  modified  ISDC  procedure  can  be  summarized  as  follows. 

1.  Use  the  ISDC  procedure  of  section  2.3. 

2.  Modify  the  equation  of  motion  of  substructure  (k-1)  by  adding  the  equivalent  dynamic 
loading,  which  includes  the  interface  stiffness,  inertial  loading  and  interacting  force  from 
substructure  (k+1),  of  substructure  k  and  interacting  force  from  substructure  (k-2). 

3.  Design  the  controller  using  the  modified  EOM  of  substructure  (k-1). 

4.  Similarly,  modify  the  equation  of  motion  of  substructure  (k+1)  by  adding  the  equivalent 
dynamic  loading  of  substructure  k  and  interacting  force  from  substructure  (k+2). 

5.  Use  the  ISDC  procedure  of  section  2.3. 

The  simple  flowchart  for  the  general  ISDC  procedure  is  presented  in  Figure  2.2. 

2,5  ISDC  for  General  Structures 

For  more  general  structures,  any  substructure  might  connect  with  more  than  two 
substructures.  The  direct  condensation  process  will  lead  to  a  difficulty  in  formulation  of  the 
tranformation  matrix.  The  indirect  condensation  process  can  express  the  transformation  matrix 
as  in  sec.  (2.2.2)  in  which  interior  degrees  of  freedom  include  the  original  interior  d.o.f.’s  and 
all  other  boundary  d.o.f.’s  which  are  not  needed  in  the  final  expression.  For  example,  a 
substructure  k,  which  connects  m  substructures,  contains  the  original  interior  degrees  of 
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freedom  (v,)  and  the  boundary  degrees  of  freedom  (vj, )  which  connect  substructures  k  and  i 
(i=kl,k24c3,...,km).  If  the  transformation  matrix  between  and  all  other  d.o.f.’s  is  going  to 
be  found  then  the  interior  d.o.f.’s  in  Eq.(2.35)  consist  of  Va  and  Vb,  (j=k2,k3,...,km),  and  the 
boundary  d.o.f.’s  consist  of  Vki.  Similarly,  the  transformation  matrix  between  each  set  of 
boundary  d.o.f.’s  and  the  other  d.o.f.’s  can  be  found.  Hence,  the  indirect  condensation  process 
makes  the  transformation  very  simple  and  can  be  used  in  any  general  structure. 

2.6  Examples 

Two  examples,  involving  a  two-bay  truss  and  a  six-bay  truss,  are  considered  to  illustrate 
the  procedure  of  interacting  substructure  decentralized  control  outlined  in  this  chapter.  For  the 
two-bay  truss  (Fig.  2.3),  two  substructures  are  used  and  four  sensors  and  actuators  are 
colocated  at  degrees  of  freedom  1,2,5,  and  6.  The  weighting  matrices  Q  and  R  are  selected  as 
10^*1  and  I,  respectively,  where  I  is  the  identity  matrix.  Following  the  procedure  described  in 
sec.  2.3,  the  closed-loop  eigenvalues  obtained  are  shown  in  Fig.2.4  which  are  compared  with 
those  of  the  controlled  structure  designed  using  the  complete  system  model.  Most  of  the 
closed-loop  eigenvalues  are  close  for  these  two  systems.  Furthermore,  the  first  three  modes  for 
the  system  with  ISDC  design  have  better  damping  property  than  the  system  with  LQR 
controller.  The  simulation  results  for  the  dynamic  responses  of  degrees  of  freedom  1,2,5,  and  6 
under  a  unit  impulse  at  degrees  of  freedom  1  and  2  is  shown  in  Figs.2.5  to  2.8.  These  results 
indicate  that  the  interacting  substructure  decentralized  control  algorithm  leads  to  a  better 
dynamic  response.  It  doesn’t  mean  that  the  ISDC  approach  is  better  than  the  LQR  design 
because  the  performance  index  might  not  be  better  for  the  system  with  the  ISDC  design. 

For  the  six-bay  truss  (Fig.2.9),  three  substructures  are  used  and  six  sensors  and  actuators 
are  colocated  at  degrees  of  freedom  1,2,5,6,17,  and  18.  The  weighting  matrices  Q  and  R  are 
chosen  as  10^*1  and  I,  respectively.  The  global  controller  is  designed  using  Eqs.(2.37)-(2.66). 
The  closed-loop  eigenvalues  are  plotted  in  Fig.2.10  and  compared  with  those  of  the  LQR 
controller  designed  using  the  complete  system  model.  The  ISDC  design  gives  a  more 
conservative  result  than  the  LQR  controller  with  complete  system  model  for  lower  frequency 
modes.  Figures  2.1 1  to  2.16  show  the  simulation  results  of  degrees  of  freedom  1,2,5,6,17,  and 
18  under  a  unit  impulse  at  degrees  of  freedom  1  and  2  (horizontal  and  vertical  directions  at  tip 
of  truss).  A  better  dynamic  response  is  shown  for  the  controller  designed  by  the  ISDC 
approach.  The  input  forces  of  the  six  actuators  are  illustrated  in  Figures  2.17  to  2.22.  The 
maximum  amplitudes  of  the  first  four  actuators  are  almost  the  same,  while  those  of  the  fifth  and 
the  sixth  actuators  of  the  ISDC  controller  are  about  15%  to  25%  higher  than  those  of  LQR 
controller  designed  with  complete  system  model. 

t 

The  performance  index,  J(x^Qx  +  u^Ru)dt,  is  presented  in  Fig. 2.23  for  both  the  systems. 

0 

The  performance  index  of  the  system  with  ISDC  controller  is  higher  than  that  of  the  system 
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with  LQR  controller  by  17.7%  for  the  case  of  impulse  input.  The  second  case  of  simulation 
involves  applying  the  initial  displacements  to  the  various  degrees  of  freedom,  which  are 
obtained  by  applying  a  static  force  at  node  I  in  vertical  direction  such  that  the  displacement  at 
node  1  is  unity.  The  dynamic  responses  of  degrees  of  freedom  1,2,5,6,17,  and  18  are  shown  in 
Figures  2.24  to  2.29.  The  dynamic  response  of  the  system  with  ISDC  controller  is  a  little  better 
than  that  of  the  system  with  LQR  controller.  Figures  2.30  to  2.35  show  the  input  time  history 
of  both  systems.  They  show  that  the  ISDC  approach  induces  smaller  input  forces  for  the  first 
three  actuators  but  higher  input  forces  for  the  last  three  actuators.  In  Fig.  2.36,  the  performance 
index  of  the  system  with  ISDC  controller  is  about  20.1%  higher  than  that  of  the  sy.stem  with 
LQR  controller. 

For  the  six-bay  truss,  the  actuators  are  rearranged  at  d.o.f.’s  1,2,5,6,21,  and  22  and  three 
substructures  are  arranged  as  shown  in  Fig.  2.9.  Here  the  second  substructure  contains  no 
actuator,  and  the  design  procedure  of  sec.  2.4  is  then  used  to  solve  the  control  design  problem. 
The  resulting  closed-loop  eigenvalues  are  compared  with  those  of  the  optimal  design  in  Figure 
2.37. 

The  ISDC  procedures  with  direct  and  indirect  condensation  approaches  are  compared  in 
Figure  2.38  for  the  six-bay  truss  described  in  Fig.  2.9.  The  closed-loop  system  eigenvalues 
obtained  using  the  indirect  condensation  approach  are  closer  to  those  of  the  optimal  design  than 
those  using  the  direct  condensation  approach.  This  shows  that  the  indirect  condensation 
approach  is  better  than  the  direct  condensation  approach. 

2.7  Concluding  Remarks 

The  interacting  substructure  decentralized  control  design  procedure  has  been  developed  for 
large  flexible  structural  control.  This  method  is  based  on  the  physical  coordinates  instead  of 
modal  coordinates.  Thus,  the  natural  frequency  and  modal  matrix  information  is  not  required 
and  no  model  data  error  is  introduced  into  the  system.  By  using  the  concept  of  interacting 
substructures,  the  size  of  the  mathematical  model  can  be  reduced  without  using  model 
reduction  techniques  and  hence  the  "spillover  problem"  can  be  avoided.  During  the  nodal 
condensation  process,  two  approaches,  direct  and  indirect  condensation  processes,  are 
considered.  Although  no  modal  data  is  involved  in  the  direct  condensation  approach,  it  can 
only  be  applied  to  serial  structures.  While  the  indirect  condensation  involves  substructure 
modal  data,  it  can  be  applied  to  any  general  structure.  As  can  be  expected,  the  design  of  the 
ISDC  using  indirect  condensation  approach  will  be  closer  to  the  optimal  design  by  increasing 
the  number  of  modes  used  in  finding  the  transformation  matrix,  which  can  be  done  by 
increasing  the  number  of  boundary  degrees  of  freedom  when  the  substructures  are  defined. 
Generally  speaking,  the  computational  efficiency  of  the  ISDC  approach  will  be  more 
significant  when  it  is  applied  to  very  large  flexible  structures. 
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Figure  2.1  Substructures  in  Series 


Figure  2.2  Genersl  ISDC  Procedure  for  Serial  Large  Flexibie  Structures 
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Note:  Numbers  represent  degree  of  freedom 

Numbers  wiA  circle  represent  the  node  numbers 


Figure  2.3  Two-Bsf  Tni» 
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Legend  : 

X  :  Eif enraluee  of  open*loop  structure 
A :  Eigenmiuei  of  eloeed-loop  system  with 
controller  designed  with  complete  model 
O  •  ElgenTiles  of  elosed«loop  system  with 

controller  designed  with  substructure  model 


Hgure  2.4  Eigenvalues  of  Two-Bay  Truss 
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:  Controller  designed  with  substructure  model 
:  Controller  designed  with  complete  model 


Figure  2.5  Displacement  of  2-ba7  Truss  at  d.o.f.  1  (impulse  at  tip) 
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■  —  :  Controller  designed  with  substructure  model 
. :  Controller  designed  with  complete  model 

Figure  2.8  Displacement  of  2-ba7  Truss  at  d.o.f.  2  (impulse  at  tip) 
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:  Controller  designed  with  substructure  model 
:  Controller  designed  with  complete  model 


Figure  2.7  Displacement  of  2-bs7  Truss  at  d.o.f.  5  (impulse  at  tip) 
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:  Controller  designed  with  substructure  model 
:  Controller  designed  with  complete  model 


Figure  2.8  Displacement  of  2-ba7  Truss  at  d.o.f.  6  (impulse  at  tip) 
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Figure  2.8  Six-Bay  Tru« 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.11  Displacement  of  9-ba]r  Truss  at  d.o.f.  1  (impulse  at  tip) 


:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.12  Displacement  of  O-bay  Truss  at  d.o.f.  2  (impulse  at  tip) 


:  Coatrotler  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.13  Displacement  of  6-bay  Truss  at  d.o.f.  5  (impulse  at  tip) 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.14  Displacement  of  6-ba7  Truss  at  d.o.f.  6  (impulse  at  tip) 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  vrith  complete  model 


Figure  2.15  Displacement  of  6>bay  Truss  at  d.o.f.  17  (impulse  at  tip) 


:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.16  Displacement  of  6-bay  Truss  at  d.o.f.  18  (impulse  at  tip) 


:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.19  Input  Force  of  Actuator  3  (d.o.f.  5)  for  8-bay  Truss  (impulse) 


48 


0.00  0.50  1.00  1.50  2.00  2.50  3.00  3.50  4.00 

Time  (sec) 


:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.20  Input  Force  of  Actuator  4  (d.o.f.  8)  for  6-bay  Truss  (impulse) 
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Ficure  2.21  Input  Force  of  Actuator  5  (d.o.f.  17)  for  6-bay  Truss  (impulse) 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.25  Displacement  of  O-bay  Truss  at  d.o.f.  2 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.20  Displacement  of  0-bay  Truss  at  d.o.f.  5 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figaro  2.27  Displacement  of  0-bay  Trass  at  d.o.f.  6 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.28  Displacement  of  0-bay  Truss  at  d.o.f.  17 
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:  ControUer  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.29  Displacement  of  d-bay  Tnua  at  d.o.f.  18 
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Figure  2.30  Input  r  orce  of 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.31  Input  Force  of  Actuator  2  (d.o.f.  2)  for  S-bay  Truss 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.32  Input  Force  of  Actuator  3  (d.o.f.  5)  for  6-bay  Truss 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.33 


I 


Input  Force  of  Actuator  4  (d.o.f.  8)  for  O-bay  Tiuss 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.34  Input  Force  of  Actuator  5  (d.o.f.  17)  for  6-bay  Truss 
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:  Controller  designed  using  ISDC  approach 
:  LQR  Controller  designed  with  complete  model 


Figure  2.36  Performance  Index  ts  Time  for  6>ba7  Truss 
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Figure  2.37  Eigenyelues  of  O-bmy  Truss  (first  12  pmrs) 
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Figure  2.38  Eigenvalues  of  6-bay  Truss  with  Different  Controllers  (first  12 


CHAPTERS 


DUAL  PASSIVE/ACTIVE  CONTROL 


3.1  Unified  Passive  Damper  Design 

Usually,  discrete  passive  dampers  are  added  to  flexible  structures  at  specific  locations  to 
increase  the  stability  of  the  whole  structure  as  well  as  to  improve  the  robusmess  of  the  closed- 
loop  system.  In  this  section,  an  approach  is  developed  for  passive  danqier  design.  The  basic 
idea  is  to  design,  first  a  feedback  controller  for  the  original  structure,  and  then  a  set  of  passive 
dampers  to  approximate  the  behavior  of  the  feedback  controlla*.  The  purpose  of  this  passive 
damper  design  procedure  is  to  use  the  well  developed  control  design  techniques  for  designing  a 
near-optimal  set  of  dampers  in  the  sense  of  minimizing  the  performance  index  of  the  closed- 
loop  system.  The  performance  index  is  defined  as  the  sum  of  the  vibration  energy  and  the  input 
energy  of  actuators.  The  input  energy  of  the  actuators  can  be  considered  to  be  proportional  to 
the  weight  of  the  damping  devices  in  the  unified  passive  damper  design. 

For  flexible  structures,  the  equations  of  motion  can  be  written  in  the  form  of  Eqs.(2.1)  and 
(2.2).  If  the  optimal  linear  quadratic  regulator  is  used  for  the  controller  design,  the  control  gain 
can  be  obtained  by  solving  the  algebraic  Riccati  equation  (2.6)  and  Eq.(2.S).  Partitioning  the 
Riccati  solution  matrix  P  into  four  nxn  submatrices  (where  n  denotes  the  dimension  of  the 
displacement  vector  v)  as 


P= 

Pl2  P22 

and  substimting  Eqs.(2.5)  and  (3.1)  into  Eq.(2.2)  leads  to  the  closed-loop  system 

x  =  Acix 


where 


(3.1) 


(3.2) 


Aci  = 


0  I 

-NT  1  DR"*  D'*'M-'*'P,2-M-*  K  -NT*  DR"*  D'*'M-'*'P22-M- *  C 


(3.3) 


Equation  (3.2)  can  be  rewritten  as 
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Mx  +  (C+AQx  +  (K+AK)x  =  0 

(3.4) 

where 

« 

ak=dr-‘d'*'m-'^Pi2 

(3.5) 

and 

ac=dr-*d'*‘m-'^P22 

(3.6) 

The  changes  in  C  and  K,  namely,  AC  and  AK,  brought  by  the  active  controller  can  also  be 
partly  achieved  using  a  set  of  discrete  dampers  and  springs.  When  m  dampers  (Ci)  and  p  spring 
(kj)  at  several  specified  locations  are  used,  the  damping  and  stiffness  matrices  of  the  system  are 
modified  by  AC  and  AK  respectively.  The  values  of  Ci  and  k;  can  be  determined  by  using  the 
least  squares  method  such  that  the  differences  between  AC  (AK)  and  AC  (AK)  are  minimized. 
The  least  squares  problem  is  formulated  and  solved  as  follows: 

1.  Write  AC,  AC,  AK  and  AK  in  vector  form  Cl,  Cl,  K1  and  K1  (e.g.  Clj  =  AQj  and  l=(i- 

l)*n+U  _ 

2.  Write  Cl  and  K1  into  a  constant  matrix  multiplying  an  unknown  vector  such  that 


Cl 


C1  =  A 


C2 


Cm 


(3.7) 


K1  =  B 


kz 

kp 


and  A,  B  are  n^*m  and  n^*p  coefficient  matrices. 

3. 

Min  |ia-Cl|| 

Cl  '"C. 


(3.8) 


(3.9) 


Min  IIKl- Kill  (3.10) 

k,"kp 
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Since  AC  and  AC  are  symmetric  matrices,  we  can  reduce  the  order  of  the  problem  in  Eq. 
(3.9)to  m*n(n+l)/2,  as 


Min  IJA' 

‘  Cm 


CiC2>..C|]} 


(3.11) 


where  Cl'  is  the  upper  triangular  matrix  of  AC  in  vector  form,  and  A'-|ciC2  •  ■  c^  i-  is 


4. 


the  upper  triangular  matrix  of  AC  in  vector  form.  Since  AK  is  a  nonsymmetric  matrix, 
£q.(3.10)  cannot  be  reduced. 

The  Moore-Penrose  inverse  [14]  is  used  to  solve  Eqs.(3. 1 1)  and  (3. 10)  as 

r  'iT 


c,C2...c™f  =A'^C1' 


(3.12) 


kik2...kp 


T 

=B‘"K1 


(3.13) 


where 


A'+ = vr-u'^ 


(3.14) 


Z*"  =diag[ai’*  ai*  *  •  *  <J7^0  •  0] 


A'  =  UZV^  (by  singular  value  decomposition) 


B"'  =  ViZtuT 


(3.15) 


B  =  Ui  Zi  v7  (by  singular  value  decomposition) 

5.  In  this  step,  some  adjustments  have  to  be  made.  Usually,  kj  and  cj  obtained  in  step  (4) 
are  not  all  positive  semidefinite  (i.e.,  not  all  realizable);  the  negative  elements  should  be 
set  equal  to  zero.  If  the  values  of  springs  are  very  very  small  compared  to  the  structural 
stiffness,  the  designer  can  neglect  them  without  much  effect. 
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To  get  a  more  general  approximation,  the  constrainted  optimization  techniques  can  be 
applied  in  the  presence  of  constraints  on  design  variables  (e.g.,  size  and  weight  constraints  on 
dampers  and  stiffnesses).  For  unconstrained  approximation,  the  least  square  approach  gives 
the  same  result  as  that  of  a  nonlinear  optimization  approach  but  the  former  will  save  much 
computadonal  time. 

3.2  Dual  Passive/ Active  Control 

As  we  know,  passive  control  is  unconditionally  stable  for  a  stable  open  loop  system  and  is 
robust  and  reliable.  Unfortunately,  passive  control  fails  to  meet  stringent  performance 
requirements  which  are  required  in  most  flexible  space  structures.  On  the  other  hand,  active 
control  meets  the  stringent  performance  requirements  if  there  are  no  energy  and  power 
limitations  and  if  a  perfect  system  model  is  available.  For  practical  structures,  the  model 
reduction  procedures  are  to  be  used;  also  there  will  be  limitations  on  energy  and  power.  In 
extreme  cases,  the  spillover  problem  caused  by  mode  truncation  will  induce  an  unstable 
closed- loop  system.  Thus,  dual  passive/active  control  approach  can  be  used  as  a  compromise 
approach  to  meet  the  performance  requirements  as  well  as  to  make  it  more  tolerant  to 
uncenainties  in  system  design  and  disturbances.  Also,  the  dual  passive/active  control  design 
leads  to  a  more  reliable  control  with  a  smaller  production  cost. 

Gehling,  1986  shows  that  if  passive  damping  is  designed  into  the  system,  the  amount  of 
damping  achievable  in  a  flexible  space  structure  will  reduce  the  requirements  of  active  control 
design.  Savings  can  be  realized  in  terms  of  the  number  of  control  system  components  and 
energy  expenditure  for  vibration  control. 

The  procedure  of  passive  damper  design  (sec.  3.1)  coupled  with  the  interacting  substructure 
decentralized  control  (Chap.  2)  can  be  described  briefly  as  follows.  The  first  step  in  dual 
passive/active  control  design  is  to  use  a  well  developed  control  technique  to  design  an  active 
controller  and  use  an  approximation  process  to  design  a  set  of  preassigned  dampers  as 
described  in  Sec.  3.1  The  second  step  is  to  modify  the  original  structure  with  the  additional 
passive  dampers  and  rewrite  the  system  equations  of  motion.  The  sensor/actuator  location 
should  be  preassigned  and  then  the  state  equation  can  be  wrinen.  Finally  the  interacting 
substructure  decentralized  control  approach  is  applied  to  design  the  active  control  system. 

3.3  Examples 

Four  examples  are  considered  to  illustrate  the  procedures  outlined  in  this  chapter.  The  first 
one  deals  with  the  passive  damper  design  using  the  unified  passive  damping  design  (UPD) 
method  on  a  two-bay  truss.  The  second  and  third  examples  deal  with  the  active  control  using 
the  ISDC  method  on  two-bay  and  six-bay  trusses,  respectively.  The  fourth  example  deals  with 
the  multilevel  passive/active  control  using  the  UPD  and  the  ISDC  methods  on  a  two-bay  truss. 
For  original  stniCtures,  no  damping  is  assumed  and  hence  the  open-loop  eigenvalues  of  the 
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structures  are  all  located  on  the  imaginary  axis. 

3.3.1  Passive  Damper  Design 

A  two-bay  truss  (Fig.  2.3),  with  prespecified  arrangement  of  damping  devices,  is  used  to 
illustrate  the  passive  damper  design  methodology.  Two  types  of  damping  arrangements,  as 
shown  in  Fig.  3.1,  are  considered  in  this  example.  The  weighting  matrices  Q  and  R  are  chosen 
as  10^*1  and  I  (I  is  the  identity  matrix),  respectively.  Following  the  UPD  procedure,  the 
damping  coefficients  have  been  determined  as  indicated  in  Table  3.1.  The  eigenvalue  pairs  of 
the  modified  structures  (the  original  structiire  with  passive  dampers)  are  compared  with  those 
of  the  closed-loop  controlled  structure  with  active  controller  designed  by  the  LQR  procedure. 
The  first  configuration  gives  an  acceptable  approximation  to  the  closed-loop  system  with  active 
control  except  for  the  first  three  eigenvalue  pairs  (see  Fig.  3.2).  The  first  four  eigenvidue  pairs 
of  the  second  configuration  are  not  very  close  to  those  of  the  closed-loop  system  with  active 
control  (Fig.  3.2).  This  shows  that  the  passive  damper  design  procedure,  described  in  this 
chapter,  can  improve  the  system  dynamic  response  very  wed  in  higher  frequency  modes. 
However,  it  is  impossible  to  obtain  a  passive  damper  design  which  gives  the  exact  dynamic 
response  of  the  actively  controlled  structure.  Hence,  the  combined  passive  and  active  control 
system  is  necessary  for  practical  structures  with  stringent  dynamic  requirements.  The  damping 
coefficients  vary  with  the  weighting  matrices  chosen.  The  relationship,  for  the  second 
configuration  (Fig.  3.1(b))  of  damping  arrangement,  is  presented  in  Table  3.2.  The  input 
weighting  matrix  R  is  the  identity  matrix  and  the  output  weighting  matrix  Q  is  the  identity 
matrix  multiplied  by  a  scalar  coefficient  q.  Tables  3.2  and  3.3  show  that  the  change  in  damping 
coefficients,  as  well  as  the  change  in  the  real  part  of  eigenvalues  of  the  modified  structures,  is 
approximately  proportional  to  the  square  root  of  the  change  in  the  scalar  coefficient  q. 

3  J.2  Dual  Passive/ Active  Control 

The  unified  passive  damping  design  procedure  (Fig.  3.2)  gives  a  good  approximation  to  the 
LQR  actively  controlled  structure  except  for  the  first  few  modes  and  the  interacting 
substructure  decentralized  control  (ISDC)  method  gives  a  better  compensation  for  the  first  few 
ryodes.  Hence,  a  dual  passive/active  control  design  for  LFS  with  the  first-level  control 
designed  by  the  unified  passive  damping  approach,  and  the  second-level  control  designed  by 
uping  ISDC,  is  expected  to  be  more  promising  in  obtaining  a  better  performance.  The  two-bay 
truss  shown  in  Fig.  2.3  is  considered  to  demonstrate  the  multilevel  passive/activc  control 
design  procedure.  The  results  of  the  two- bay  truss  configuration  shown  in  Fig.  3.1(b),  obtained 
v^th  only  the  passive  damper  design  and  weighting  matrix  Q=  10^*1,  are  used  as  the  first  level 
passive  control  design  in  this  example.  A  second  level  active  control  design  is  obtained  by 
using  the  ISDC  procedure  with  weighting  mattices  Q=10^*I  and  R=I  based  on  the  passively 
camped  structure.  Two  substructures  are  used  and  four  sensors  and  actuators  are  colocated  at 
degrees  of  freedom  1,2,5,  and  6.  The  eigenvalue  pairs  of  (i)  the  uncontrolled  structure,  (ii)  the 
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controlled  structure  with  the  controller  designed  with  the  complete  model,  (iii)  the  controlled 
structure  with  controller  designed  by  the  ISDC  method,  and  (iv)  the  structure  designed  with  the 
dual  passive/active  control ,  are  shown  in  Fig.  3.3. 

Most  of  the  eigenvalues  of  controlled  structure  with  UPD/ISDC  controller  move  to  the  left 
compared  to  those  of  the  closed-loop  system  with  ISDC  controller.  This  indicates  that  the 
stability  margin  of  the  closed-loop  system,  for  the  higher  frequency  modes,  is  increased  by 
using  the  dual  passive/active  control  strategy.  Although  the  change  is  small,  the  reliability  and 
stability  robusmess  of  the  system  will  be  significantly  improved.  The  relationship  between  the 
change  in  the  eigenvalues  of  the  dual  UPD/ISDC  controlled  structure  and  that  in  the  output 
weighting  matrix  is  also  illustrated  in  Table  3.4. 

The  simulation  results  of  the  two-bay  truss  with  unit  impulse  at  tip  (in  both  directions  at 
node  1)  for  different  controllers  are  presented  in  Figs.  3.4  to  3.11.  The  displacement  outputs 
(Figs.  3.4  to  3.7)  from  the  sensors  appear  to  be  identical  for  systems  with  ISDC  controller  and 
dual  UPD/ISDC  controller.  Figures  3.8  to  3.11  show  that  the  input  forces  of  ISDC  controller 
are  larger  than  those  of  dual  UPD/ISDC  controller.  The  time  variation  of  the  performance 
index  is  presented  in  Fig.  3.12.  The  performance  index  of  the  dual  UPDASDC  controller  is 
reduced  by  10.7%  compared  to  that  of  LQR  controller  and  by  28.14%  compared  to  that  of 
ISDC  controller.  This  indicates  that  a  great  improvement  in  performance  has  been  achieved  by 
the  dual  passive/active  control  design. 

3.4  Concluding  Remarks 

A  unified  passive  damper  design  procedure  is  developed  based  on  the  well  developed 
control  design  techniques  and  approximation  techniques.  The  design  variables  are  the 
weighting  matrices  R  and  Q  for  LXJR  design  or  control  gain  for  other  control  design 
techniques.  Appropriate  location  selection  of  damping  devices  is  very  important  for  the  unified 
passive  damper  design  procedure.  For  the  case  when  the  weighting  matrices,  R  and  Q,  are 
assumed  to  be  the  identity  matrix  and  the  identity  matrix  multiplied  by  a  scalar  q,  the  change  in 
the  damping  coefficients,  for  a  preassigned  set  of  damping  devices,  varies  as  the  square  root  of 
the  change  in  the  scalar  coefficient  q.  The  design  of  constrained  viscoelastic  layer  treatment,  a 
kind  of  continuous  damper,  can  also  be  achieved  by  this  method  using  optimization  techniques. 
The  dual  passive/active  control  strategy,  using  unified  passive  damper  design  and  interacting 
substructure  decentralized  control  approaches,  is  applied  for  a  two-bay  truss.  The  result  shows 
that  the  passive  control  significantly  improves  the  reliability  of  the  closed-loop  system  and 
reduces  the  active  control  energy  consumption  under  the  same  performance  requirements. 
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Table  3.1  Damping  Coefficients  of  the  UPD  for  Two^Bay  Truss  (Qs=  10**1) 


Dsmoinc  Coeff.  (lb*s 

/in^ 

I 

n 

Cl 

31.6271 

15.8116 

Ct 

21.0015 

21.0015 

C3 

31.8271 

15.8116 

C4 

31.6267 

15.8151 

cs 

21.0871 

10.5457 

Cfl 

31.6267 

15.8151 

Table  3.2  Damping  Coefficients  of  the  UPD  for  Two-Bay  Truss  in 
Configuration  II 


Damnins  Coeff.  for  Configuration  II 

Matrix  Q 

10*1 

10**1 

10**1 

10^*1 

10**1 

Cl 

1.5812 

15.8116 

50.0006 

158.1160 

cj 

6.6607 

21.0015 

66.6072 

210.9151 

C3 

1.5812 

15.8116 

50.0006 

158.1160 

C4 

1.5815 

KImT  m 

15.8151 

50.0117 

158.1510 

Cs 

3.3348 

10.5457 

33.3485 

105.4572 

_ Cft. 

5.0012 

15.8151 

50.0117 

158.1510 
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Table  3.3  Eigenyaluea  of  Modified  Two^Bay  Truia  with  Damping  Coefficients 
of  Table  3.2 


E 

leenvalues  of  Modified  Two-Bar  Truss 

Matrix  0. 

10^*1 

10^*1 

10*  *1 

Eigenvalue 

pairs 

-0.28-l-143.97i 

-0.28-l-124.84i 

-0.23-t-123.06i 

.0.10-l-104.11i 

-0.03-f-87.27i 

-0.01-t-46.10i 

-0.02-M2.69i 

-0.002-t-14.10i 

-0.88-t-143.97i 

-0.88-f-124.84i 

-0.72-t-123.06i 

-0.32-f-104.lli 

-0.09-t-87.271 

-0.03-M6.10i 

-0.08-M2.69i 

-0.01-l-14.10i 

-2.78-t-143.93i 
-2.78-t-124.82i 
-2.26-t-123.04i 
-1.02-l-104.11i 
-0.29-t-87.27i 
-0.10-M6.10i 
-0.25-M2.69i 
-0.02-1-14.  lOi 

Table  3.4  Elgenralues  of  Cloeed-Loop  Two-Bay  Trass  with  Different 
Controllers 


Sieenvaiue  Pairs 

LQR 

ISDC 

ICSBSBm 

IcSSISiHI 

-l.S5-t-143.96i 
-1.48-t-124.73i 
-1.29-t-123.15i 
-l.02-t-104.lli 
-1.39-t-  87.26i 
-1.25-1-  46.07i 
-0.88-t-  42.68i 
•0.99-f  14.06i 

-l.73-i-144.01i 
-1.43-1-124.061 
-I.10-fl23.79i 
-1.06-fl04.17i 
-1.02 -f  87.28i 
-2.19-f  46.24i 
-1.23-f  42.67i 
-1.33-f  14.12i 

-l.79-fl43.99i 
-I.3S-fl24.29i 
-1.38-fl23.57i 
.1.13-fl04.17i 
-1.05+  87.28i 
-2.14+  46.19i 
-1.27+  42.97i 
-1.34+  U.lli 

-2.06+143.971 
-1.61+124.52! 
-1.68+123.351 
-1.31+104.151 
-1.12+  87.281 
•2.06+  46.19i 
-1.36+  42.69! 
-1.35+  14.101 

-3.47+143.911 
-3.10+124.711 
-2.87+123.141 
-1.94+104.131 
-1.31+  87.281 
-1.95+  46.151 
•1.50+  42.691 
-1.37+  14.09! 
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(a) 


Damping  Arrangement  for  Two-Bay  Truss  (a)  Configuration  I,  (b) 
Configuration  II 
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X  :  Eigenvalues  of  original  structure 
A  :  Eigenvalues  of  closed-loop  system  with 
controller  designed  by  LQR 
O  :  Eigenvalues  of  modified  system  I  (Fig.  3(a)) 

□  :  Eigenvalues  of  modified  system  II  (Fig.  3(b)) 

Figure  3.2  Eigenvalues  of  Closed-Loop  Controlled  and  Modified  2-bay  Truss 


a  A 


X 


r  150.0 


Re 


Legend  : 

X  :  Eigenvalues  of  open-loop  structure 
A  :  Eigenvalues  of  closed-loop  system  with 
controller  designed  with  complete  model 
□  :  Eigenvalues  of  closed-loop  system  with 
the  ISDC  controller 

O  :  Eigenvalues  of  closed-loop  system  with 
multi-level  UPD/ISDC  controller 

Figure  3.3  Eigenvalues  of  2-bay  Truss 
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_ :  UPDASDC  two-levcl  controller 

- :  ISDC  concrol  design 

. :  LQR  active  control  with  complete  model 

Figure  3.4  Displacement  of  2-bay  Truss  at  d.o.f.  1  (impulse  at  tip) 
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_ :  UPD/ISDC  two-lcvel  controller 

- ;  ISDC  control  design 

. :  LQR  active  control  with  complete  model 

Figure  3.5  Displacement  of  2-bay  Truss  at  d.o.f.  2  (impulse  at  tip) 


_ :  UPD/ISDC  two-level  controller 

- ;  ISDC  control  design 

. :  LQR  active  control  with  complete  model 

Figure  3.6  Displacement  of  2-bay  Truss  at  d.o.f.  5  (impulse  at  tip) 
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_ :  UPDASDC  two-ievel  controller 

- :  ISDC  control  design 

. :  LQR  active  control  with  complete  model 

Figure  3.7  Displacement  of  2-bay  Trtiss  at  d.o.f.  6  (impulse  at  tip) 


0.00  0.50  1.00 


1.50  2.00  2.50  3.00 
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, :  UPD/ISDC  two-levei  controller 
:  ISDC  control  design 
:  LQR  active  control  with  complete  model 


Figure  3.8  Input  Force  of  Actuator  1  of  2>ba]r  Truss  (impulse  at  tip) 


8. 


Time  (sec) 


_ :  UPD/ISPC  two-levci  conoDUer 

- :  ISDC  control  design 

. :  LQR  active  control  with  complete  model 

Figure  3.10  Input  Force  of  Actuator  3  of  2-ba7  Truss  (impulse  at  tip) 
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_ :  UPD/ISDC  two-level  controller 

- :  ISDC  control  desip 

. :  LQR  active  control  with  complete  model 


Figure  3.12  Performance  Index  of  2-bay  Truae  (impulse  at  tip) 


CHAPTER  4 


ROBUST  INTEGRATED  STRUCTURAL  AND  CONTROL  DESIGN 


4.1  Background 

There  has  been  a  dramatic  increase  in  the  past  decade  in  using  active  control  systems  to 
improve  structural  performance.  The  major  challenge  in  the  field  of  active  control  of  structures 
is  in  the  design  of  control  systems  for  very  large  space  structures.  These  structures  are  by 
nature  distributed  parameter  systems  with  multiple  inputs  (controls)  and  a  continuum  of 
outputs  (displacements).  The  finite  element  method  is  commonly  used  for  the  description  of 
these  structures.  This  is  a  source  of  parameter  errors  and  truncated  (or  reduced  order)  models  in 
the  system.  In  addition,  the  structural  propoties  of  large  space  structures  cannot  be  tested 
before  they  are  put  into  orbit  and  hence  sizeable  uncertainties  exist  in  modal  parameters. 

A  great  deal  of  research  is  currently  in  progress  on  developing  methods  for  the 
simultaneous  (integrated)  design  of  the  structure  and  the  control  system.  The  weight  of  the 
structure  was  minimized  with  constraints  on  the  distribution  of  the  eigenvalues  and/or  damping 
ratio  of  the  closed-loop  system  by  Khot  et  al.  (1985).  Miller  and  Shim  (1986)  considered  the 
simultaneous  minimization,  in  structural  and  control  variables,  of  the  sum  of  structural  weight 
and  the  infinite  horizon  linear  regulator  quadratic  control  cost.  The  sttucture/control  system 
optimization  problem  was  formulated  by  Khot,  et  al  (1986)  with  constraints  on  the  closed-loop 
eigenvalue  distribution  and  the  minimum  Frobetuous  norm  of  the  control  gains.  In  all  the 
above  works,  the  consideration  of  robusmess  of  the  control  system  has  been  ignored. 

The  parameter  variations  introduced  by  the  analysis  model,  uncertain  material  properties  or 
optimization  may  adversely  influence  the  stability  and  performance  characteristics  of  the 
control  system.  The  robustness  is  an  extremely  important  feature  of  a  feedback  control  design. 
A  robust  control  design  is  one  which  satisfactorily  meets  the  system  specifications  even  in  the 
presence  of  parameter  uncertainties  and  other  modeling  errors.  Since  the  system  specifications 
could  be  in  terms  of  stability  and/or  performance,  two  types  of  robusmess,  namely,  stability 
robusmess  and  performance  robusmess,  are  to  be  considered  in  the  design  stage. 

The  current  published  literature  on  control  system  robusmess  addresses  either  the  stability 
robustness  aspect  or  the  performance  robusmess  aspect.  Most  of  the  work  on  the  stability 
robusmess  (in  the  controls  area)  was  done  in  the  frequency  domain  using  singular  value 
decomposition  while  much  of  the  useful  research  on  performance  robusmess  was  carried  out  in 
time  domain  using  sensitivity  approaches.  Design  smdies  that  treated  the  stability  robusmess 
aspect  in  time  domain  and  studies  which  combined  both  stability  robusmess  and  performance 
robusmess  in  the  design  process  have  been  scarce.  The  recent  developments  in  the  area  of 
robust  multivariable  control  theory  have  been  summarized  by  Ridgely  and  Banda  (1986).  The 
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stability  robustness  of  linear  systems  was  analyzed  in  the  time  domain  in  Yedavalli  et  al.(198S) 
wherein  a  bound  on  the  perturbation  of  an  asymptotically  stable  linear  system  was  obtained  to 
maintain  stability  using  Lyapunov  matrix  equation  solution.  In  Kosut  et  al.(1983),  singular 
value  robusmess  measures  were  used  to  compare  the  performance  and  stability  robusmess 
properties  of  different  control  design  techniques  in  the  presence  of  residual  modal  intmction 
for  a  flexible  spacecraft  system.  The  importance  of  robusmess  considerations  in  the  design  of 
flexible  space  structures  was  discussed  by  Hoehne  (1985).  Gordon  and  Collins  (1985) 
presented  a  direct  design  method  for  solving  the  problem  of  robusmess  to  cross-coupling 
perturbations  by  treating  the  feedback  gains  as  design  variables.  Their  method  makes  use  of 
nonlinear  programming  techniques  along  with  a  time  domain  pole  placement  procedure.  A 
technique  for  the  improvement  of  stability  robusmess  by  shaping  the  singular  value  spectrum 
using  constrained  optimization  methods  was  described  in  Mukhopadhyay  (1985). 

4.2  Robustness  Analysis  in  Controlled  Structures 

Stability  robustness  and  performance  robustness  are  two  very  important  properties  which 
give  the  ability  of  the  closed-loop  system  to  maintain  its  asymptotic  stability  and  performance 
requirements  under  penurbations  or  errors  in  the  system  model  or  controller  parameters. 
Hence,  the  measures  of  the  stability  and  performance  robusmess  are  important  The  stability 
robusmess  index,  a  measure  of  closed-loop  system  stability  under  uncertainty  or  modeling 
errors,  is  defined  as  the  summation  of  weighted  variations  of  system  eigenvalues  and  can  be 
written  as 

P«  =  Zwi|Xi-Xil  (4.1) 

1=1 


where  is  the  stability  robusmess  index,  and  and  Xj  are  the  closed-loop  eigenvalues  of 
nominal  system  and  perturbed  system,  respectively.  The  weighting  coefficient  Wi,  denoting  the 
importance  of  eigenvalue  pair  i,  is  determined  by  the  horizontal  distance  of  eigenvalue  to 
imaginary  axis  in  complex  plane: 


1 

S*real(Xj) 


(4.2) 


in  which 


o  _  ^  1 

Jreal(Xi) 


(4.3) 


In  addition  to  stability  robusmess,  it  is  desirable  to  retain  the  performance  unchanged  when  the 
design  variables  arc  changed.  Since  the  performance  cannot  remain  the  same  ,  a  performance 


89 


robustness  index  is  defined  to  present  the  performance  robustness.  From  optimal  linear 
quadratic  regulator,  if  P  is  the  steady-state  solution  of  matrix  Riccati  equation,  the  performance 
index  is  given  as  a  function  of  initial  state  vector  as 

J  =  xJPxo  (4.4) 


The  performance  index  for  perturbed  system  with  the  same  controller  can  be  written  as 

T=  xj^  (4.5) 


where  P  satisfies  the  Lyapunov  equation 

AgjP  +  PAci  +  G^RG  +  Q  =  0 


(4.6> 


in  which,  A^  is  perturbed  closed-loop  system  matrix: 

Aci  =  A-»-BG  (4.7) 

R  and  Q  are  positive  definite  input  and  positive  semi-definite  output  matrices  and  G  is  the  gain 
matrix  for  controller: 

G  =  >-R-^B’^P  (4.8) 


and  P  satisfies  the  relation 


The  performance  robusmess 
index: 


A'  P  -I-  PA  -  PBR"^B’*'p  -H  Q  =  0 


index  is  defined  as  the  error  ratio  of 


J-J 


X?(P-P)Xo 

xJPxo 


(4.9) 

steady-state  performance 

(4.10) 


It  has  been  pointed  out  by  Kleinman  and  Athans  (1968)  that  the  optimal  solution  of  Eq.(4.4)  in 
general  depends  on  the  initial  state  xq.  This  result  is  not  very  useful  since  the  initial  state  is  not 
always  known.  The  effect  of  xq  is  averaged  out  by  assuming  that  xo  is  a  uniformly  distributed 
random  vector  whose  covariance  is  given  by  the  identity  matrix.  The  trace  of  P  is  proportional 
to  the  expected  value  of  J.  Hence  Eq.(4.10)  can  be  written  as 
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(4.11) 


Pi*  ~ 

For  both  nominal  and  perturbed  systems,  good  dynamic  response  can  be  achieved  if  the  real 
part  of  every  eigenvalue  is  restricted  to  be  smaller  than  a  specified  value  as  Re(A,)  <  -a,  with 
a>0.  The  LQ  regulator  can  be  modified  (Jacobson  et  al.,  1980)  as  follows  with  the  requirement 


of  (A,B)  being  controllable  and  (A,C)  being  observable. 

X  =  (A  +  al  )x  +  Bu  =  (A+aI-BG)x  (4. 12) 

u  =  -R"‘  B'^K.x  =  -Gx  (4.13) 

with  Ka  satisfying  the  equation 

(A+aD'^K,  +  Ka(A+aI)  -  K,BR-^  B’^K,  +  C’^QC  =  0  (4. 14) 

The  characteristic  equation  of  Eq.(4.12)  is  given  by 

det  [(A  +  al  -  BG)  -  AJ)=det  [A  -  BG  -  X'l]  =  0  (4. 15) 


with  X  =  X,  -  a.  Since  A  is  assumed  to  be  stable,  then  Re(X )  <0,  and  hence  Re(X,)  <-a.  Thus, 
the  system  with  controller,  Eq.(4.13),  provides  a  guaranteed  stability  margin. 

4  J  Multiobjective  Design  Problem 

The  stability  robusmess,  the  performance  robusmess  and  the  total  stmctural  weight  are 
considered  as  the  objective  functions  in  this  work.  The  cross  sectional  areas  of  the  members 
are  treated  as  the  design  variables.  The  first  objective  function,  the  stability  robustness  index 
(Psf),  describes  the  relative  stability  of  the  system  when  the  design  variables  change  by  a 
specified  amount.  It  is  assumed  that  the  controller  gains  are  such  that  the  condition  for  the 
stability  of  the  system  is  satisfied  and  thus  the  closed-loop  system  matrix  of  the  perturbed 
system  is  still  stable.  According  to  this  definition,  P5r=0  corresponds  to  a  highly  robust  system 
from  the  stability  point  of  view.  However,  P„  will  not  attain  the  value  zero  because  of  the 
presence  of  perturbations  in  the  design  variables.  The  second  objective  function,  the 
performance  robusmess  index  of  the  system  (Ppr),  is  defined  by  Eq.(4.11).  Here  also,  |3pr=0 
corresponds  to  a  highly  robust  system  from  the  performance  point  of  view.  The  value  ^=0 
will  not  be  attained  in  practice  because  of  the  presence  of  perturbations  in  the  design  variables. 
The  third  design  objective  function,  the  total  structural  weight,  is  given  by 

f3(x)  =  ZPiliAi  (4.16) 

i=l 

with  Pi,  Ij  and  Ai  denoting  the  density,  length  and  cross-sectional  area  of  the  i*  member 
respectively,  and  N  representing  the  number  of  members  in  the  truss  structure.  The  following 
constraints  are  used  during  the  optimization  proceduie; 


Tr(P)-Tr(P) 

Tr(P) 
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1.  Upper  and  lower  bounds  on  the  design  variables. 

2.  Stability  requirement,  i.e.,  the  requirement  of  the  real  parts  of  the  eigenvalues  of  the 
closed-loop  system  to  be  negative. 

3.  Lower  and/or  upper  bounds  on  the  natural  frequencies  of  vibration  of  the  structure. 

In  some  cases,  the  closed-loop  damping  ratios  of  the  system  may  have  to  be  constrained; 
however,  these  are  not  considered  in  this  work. 

4.4  Multiobjective  Optimization  Techniques 


The  three  objective  optimization  problem  formulated  is  solved  using  the  utility  function, 
the  lexicographic  and  the  goal  programming  methods.  The  utility  function  method  involves  the 
solution  of  the  following  problem  (Rao,  1984A): 


Min  U(x)  =  Y,  Wifi(x) 


(4.17) 


subject  to 


gj(x)^0,  j=l,2,...,m 


h  ^ 

where  w^  is  the  weight  of  the  i*^  objective  function  fj  and  Usually,  the  scales  and  units 

i=l 

of  different  objective  functions  are  different.  Hence  a  suitable  normalization  process  has  to  be 
used  in  constructing  the  objective  functions  of  Eq.(4.17).  A  convenient  form  is  to  define  a  new 
objective  function  F,  as 


Fi(x)  = 


fi(x)-fr(xr) 

fi(Xi) 


(4.18) 


where  x*  is  the  optimal  solution  for  individual  objective  function  fj,  and  f*(x*)  is  optimal 
function  value  of  objective  function  i.  Redefine  the  objective  function  of  Eq.(4.17)  as 


Min  U(x)= 

i=l 


(4.19) 


In  the  lexicographic  method,  the  objectives  are  ranked  in  order  of  importance  by  the 
designer.  The  preferred  solution  obtained  by  this  method  is  one  which  minimizes  the  objectives 
starting  with  the  most  important  one  and  proceeding  according  to  the  order  of  importance  of 
the  objectives.  Let  the  subscripts  of  the  objectives  indicate  not  only  the  objective  function 
number,  but  also  the  priorities  of  the  objectives.  Thus  Fi(x)  and  F|((x)  denote  the  most  and  the 
least  important  objective  functions,  respectively.  Then  the  first  problem  is  formulated  as 


Min  Fi(x) 


(4.20) 


subject  to 
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gj(x)^0,  j=l,2,...,m 

and  its  solution  x *  and  F*  =  Fj  (xj )  arc  obtained.  Then  the  second  problem  is  formulated  as 

Min  F2(x)  (4.21) 

subject  to 

gj(x)  0  ,  j=l,2,...,m  .and 


gm+i  (x)  =  Fi  (x)  -  Fi  (x  J ^  El 

where  Ei  is  a  small  value  compared  to  Ft  (x*).  This  procedure  is  repeated  until  all  k  objectives 
have  been  considered.  The  i*  problem  is  given  by 

Min  Fi(x)  (4.22) 

subject  to 


gj(x)  ^  0 ,  j=l,2,...,m  .and 


Sin+n(x)  ~  Fn(x)  F{i(Xn)  ^  Eo.  n— 1.2.....i  1 

The  solution  obtained  at  the  end.  i.e.  x|[  is  taken  as  the  desired  solution  x*  of  the  multiobjective 
optimization  problem. 

The  goal  programming  method  was  originally  proposed  by  Giarnes  and  Cooper  (1977)  for 
linear  optimization  problems.  The  method  requires  goals  to  be  set  for  each  objective  that  the 
designer  wishes  to  obtain.  A  preferred  solution  is  then  defined  as  the  one  which  minimizes  the 
deviations  from  the  set  goals.  Thus  a  simple  goal  programming  problem  can  be  defined  as 

Min  F(x)  =  [  X(  Fi(x)  )P  .  P  ^  1  (4.23) 

(=1 

subject  to 


gj(x)  ^  0 .  j=1.2.....m 


Fj(x)^0.  j=1.2.....k 

where  Fj(x)  =  Fj(x)-Fj(Xj* ). 

4.5  Computational  Procedure 
4.5.1  Analysis 


The  following  analysis  procedure  is  used  to  study  the  effects  of  variations  in  the  parameters 
of  the  structure  on  the  robustness  of  the  system: 
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1.  Start  with  an  initial  reference  design  of  the  structure  and  hnd  the  corresponding  plant 
matrix  A,  input  matrix  B  and  the  output  matrix  C. 

2.  Use  the  LQ  regulator  design  technique  to  find  the  optimal  control  gain  G  by  solving  the 
algebraic  Riccati  equation. 

3.  Change  the  design  parameters  by  known  percentage  values  and  find  the  corresponding  A, 
B,  and  C  matrices. 

4.  Hnd  the  stability  robusmess  index  (^•('^•1))  and  the  performance  robusmess  index 
Pp,  (Eq.(4.11)). 

5.  Repeat  steps  (3)  and  (4)  for  different  parameter  changes,  (e.g.,  nominal  design  variables, 
damping  ratio,  density  etc.) 

6.  Plot  a  graph  between  P„  or  and  the  change  in  the  parameters. 

4,5.2  Design 

The  purpose  of  design  is  to  optimize  the  actively  conoxjlled  structure  by  using  suitable 
multiobjective  optimization  techniques.  The  procedure  is  given  as  follows: 

1.  Fix)m  the  requirements  of  stress  and  deformation,  obtain  the  preliminary  design  (to  be 
used  as  the  nominal  design)  of  the  structure. 

2.  Construct  the  plant  matrix,  input  matrix,  aitd  output  matrix. 

3.  Formulate  the  multiobjective  constrained  optimization  problem. 

4.  Minimize  the  individual  objective  functions  and  find  the  respective  minima  around  the 
nominal  design. 

5.  Use  a  suitable  multiobjective  optimization  approach  to  find  a  compromise  solution. 

4.fi  Examples 
4.6.1  Two>Bar  Truss 

The  two-bar  truss  shown  in  Fig.  4.1  is  selected  for  its  simplicity.  A  nonstructural  mass  of  1 
unit  is  attached  at  node  3.  The  actuators  and  the  sensors  are  colocated  at  node  3  acting  in  x  and 
y  directions.  The  design  variables  (cross-sectional  areas  of  the  two  bars)  are  restricted  to  lie 
between  0.01  and  1.0.  The  structural  damping  ratio  is  considered  as  0.01,  Young’s  modulus  is 
assumed  to  be  10^,  and  density  is  taken  to  be  4.6.  In  the  performance  index,  the  output 
weighting  matrix  Q  is  assumed  to  be  1000.*I,  and  the  input  weighting  matrix  R  is  taken  to  be  I, 
whoe  I  is  the  identity  matrix.  The  natural  frequencies  of  the  closed-loop  controlled  structure 
are  constrained  to  lie  between  20  rad/sec  and  40  rad/sec.  For  a  stable  open  loop  system,  the 
corresponding  feedback  closed-loop  system  must  be  stable  under  the  optimum  control  law.  But 
the  stability  is  not  guaranteed  if  there  exist  disturbances  or  uncertainties  in  system  parameters. 
Hence,  additional  constraints  are  added  on  the  perturbed  closed-loop  system,  namely,  that  all 
the  eigenvalues  of  the  perturbed  closed-loop  system  are  restricted  to  have  negative  real  pans. 
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Analysis: 

Figure  4.2  shows  the  relationship  between  the  stability  robustness  index  and  the  design 
variables,  which  is  a  smooth  concave  function.  Figure  4.3  shows  the  variation  of  the 
performance  robustness  index  with  the  two  design  variables;  it  is  a  non-concave  function 
having  several  local  minima  in  the  design  space.  Figure  4.4  shows  the  variations  of  stability 
robusmess  index  and  performance  robusmess  index  with  changes  in  the  structural  dan:q>ing 
ratio  of  the  two-bar  truss.  The  stability  robusmess  index  drops  sharply  at  a  value  of  the 
damping  ratio  of  approximately  4%.  The  performance  robusmess  index  can  be  observed  to 
attain  a  minimum  value  at  a  damping  ratio  of  approximately  1%.  Figure  4.5  shows  that  the 
system  performance  index  (J)  decreases  as  the  damping  ratio  increases.  The  effect  of  the 
variations  in  Young’s  modulus  of  the  material  on  the  robusmess  indices  is  shown  in  Figure  4.6. 
The  stability  robusmess  index  monotonically  increases  as  the  Young’s  Modulus  increasing. 
The  performance  robustness  index  reduces  to  a  minimum  value  at  Young’s  modulus  (E)= 
20x10^  and  then  increases  for  larger  values  of  E.  Figure  4.7  indicates  that  J  increases  to  a 
maximum  value  at  Young’s  modulus  E=20xl0^  then  decreases  for  larger  values  of  E. 

The  relationship  of  and  Ppr  with  the  mass  density  of  the  material  is  shown  in  Fig.  4.8. 
The  stability  robusmess  index  decreases  with  an  increase  in  the  density  of  the  material.  The 
performance  robusmess  index  reduces  to  a  minimum  at  p=1.5  and  then  increases  for  higher 
values  of  p.  In  Fig.  4.9,  the  performance  index  attains  a  maximum  value  at  p=1.5  and  then 
decrease  monotonically  as  the  material  density  increases.  Figure  4.10  shows  the  variations  of 
the  stability  robusmess  index  and  the  performance  robusmess  index  with  a  change  in  the 
coefficient  of  the  output  weighting  matrix.  An  increase  in  the  coefficient  of  the  output 
weighting  matrix  implies  that  the  output  performance  is  more  important  than  the  control 
energy.  The  minimum  of  and  the  maximum  of  Ppr  occur  at  coefficient=400,  and  a  larger 
coefficient  increases  the  system  stability  robusmess  index  but  reduces  the  system  performance 
robusmess  index  when  the  coefficient  is  greater  than  400.  Figure  4.11  shows  that  the 
performance  index  reaches  the  minimum  when  the  coefficient  of  output  weighting  matrix  is 
equal  to  400  and  monotonically  increases  when  the  coefficient  is  greater  than  400. 

Design; 

The  results  of  minimization  of  the  individual  objective  functions  are  shown  in  Table  4.1. 
The  results  given  by  different  multiobjective  optimization  methods  are  shown  in  Table  4.2.  The 
first  two  columns  in  Table  4.2  correspond  to  fonnulations  #1  and  #2  of  the  utility  function 
method.  In  formulation  #1,  wj  are  set  equal  to  a  fixed  value  of  1/3  in  Eq.(2.17)  while  Wj  are 
considered  as  design  variables  in  formulation  #2.  The  last  row  of  Table  4.2  gives  the  values  of 
the  global  evaluation  function,  Fg,  which  can  be  used  as  an  index  to  compare  the  results  of 
different  multiobjective  optimization  methods.  The  global  evaluation  function  is  defined  as 
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where 


F,(X): 


(fi(x)- 0.009502) 
0.010535 


The  finite  element  model  of  the  second  example  (two-bay  truss)  is  shown  in  Fig.  4.12.  For 
this  example,  nonstructural  masses  of  magnitude  1.29  are  attached  at  nodes  1  to  4.  Each  node 
has  two  degrees  of  freedom.  The  actuators  and  sensors  are  colocated  at  nodes  1  to  4  and  are 
assumed  to  act  along  the  y-direction  only.  The  design  variables  (cross-  sectional  areas)  are 
restricted  to  lie  between  0.001  and  0.5.  The  natural  frequencies  of  the  closed-loop  system  are 
constrained  to  be  larger  than  31.62  rad/sec  (i.e.,  ca^^KXlO). 

Analysis: 

This  example  has  10  design  variables.  Since  the  display  of  functional  relations  in  10- 
dimensional  design  space  is  not  possible,  the  variation  of  the  robusmess  of  the  system  is  found 
by  uniformly  varying  the  value  of  all  the  10  design  variables.  The  results  are  shown  in  Fig. 
4.13.  This  figure  shows  the  stability  robustness  index  versus  the  value  of  the  design  variables 
when  the  permissible  change  in  the  design  vector  is  assumed  to  be  -5%.  Psr  Ppr  decrease 
slowly  with  an  increase  in  the  value  of  the  design  variables.  The  shows  a  smooth  concave 
function,  while  the  Ppr  shows  a  nonsmooth  curve.  Hence  local  minima  are  expected  in  the 
optimization  process  for  the  performance  robustness  index.  In  Fig.  4.14,  the  performance 
index  decreases  monotonically  with  increase  in  the  value  of  the  design  vector.  This  implies  that 
a  stronger  structure  will  induce  a  smaller  displacement  and  needs  lesser  control  energy  to 
obtain  good  performance. 

Design: 

The  nominal  values  of  the  design  variables  are  Xj  -0.1,  i=l  to  10.  Table  4.3  gives  the 
results  obtained  by  optimizing  the  individual  objective  functions  starting  from  the  nominal 
design.  The  results  of  different  multiobjective  optimization  methods,  namely,  the  utility 
function  method,  the  lexicographic  method  and  the  goal  programming  method  are  compared  in 


(4.25) 


Table  4.4.  The  last  row  of  Table  4.4  shows  the  global  evaluation  function,  Fg,  defined  as 

Fg(x)  =  iFi(x*) 

i=l 


where 


Fi(x)  = 


(fi(x)- 0.048694) 
0.000668 


and 


F2(x)  = 


(f2(x)- 0.008454) 
0.001588 


(f3(x)- 0.2273) 
0.02681 


4.7  Concluding  Remarks 

The  stability  robustness  index  and  the  performance  robusmess  index  defined  in  this  chapter 
are  highly  nonlinear  with  respect  to  design  variables.  The  non-concave  property  of  robustness 
indices  with  changes  in  design  variables  leads  to  difficulties  in  optimization.  As  such  one  can 
expect  to  find  only  a  local  optimum  in  the  neighborhood  of  the  starting  design  during 
optimization.  In  general,  the  local  optima  are  acceptable  since  the  starting  design  is  usually 
taken  as  the  nominal  design  which  is  expected  tt>  be  robust.  Techniques,  such  as  genetic 
algorithms  and  simulated  annealing,  might  be  alternatives  to  avoid  the  local  minima  to  obtain 
the  global  optimal  design.  The  relationships  between  the  stability/performance  robusmess 
index  and  the  various  system  parameters  have  been  determined  numerically  for  the  two-bar 
truss.  These  results  are  expected  to  be  useful  in  choosing  suitable  material  for  a  given  structure 
with  a  specified  geometry  or  weighting  coefficients  in  the  performance  index  for  controller 
design. 

A  major  advantage  of  using  nonlinear  progranuning  to  find  the  robust  control/structural 
design  is  that  it  can  be  used  with  large  permissible  changes  in  the  design  variables  and/or 
different  constraint  specifications.  Three  multiobjective  optimization  methods  have  been  used 
to  find  the  optimal  designs  of  the  illustrative  examples.  For  the  two-bar  truss,  the  utility 
function  method  with  variable  coefficients  gave  the  smallest  value  of  the  global  evaluation 
function.  For  the  two-bay  truss,  the  goal  programming  method  with  p=2  yielded  the  smallest 
value  for  the  global  evaluation  function  and  the  utility  function  method  with  variable 
coefficients  gave  the  second  smallest  value.  As  observed  in  the  investigation,  no  particular 
method  gives  the  best  solution  for  all  the  problems.  Hence,  several  methods  are  to  be  used  to 
solve  the  problem  and  find  the  best  trade-off  between  the  multiple  objectives. 
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Table  4.1 

Single  Objective  Optimization  of  Two-Bar  Truss. 

Permissible  design  variable  change  =  -5%. 

^  =  0.01,  C=103,  x(o)='qJ' 

Minimization  of 

Objective 

Psr 

Ppr 

Weight.  W 

Q 

1.0 

0. 

0. 

i=l,3 

0. 

1.0 

0. 

0. 

0. 

1.0 

0.14626 

0.15247 

0.051301 

X 

0.14626 

0.13797 

0.051301 

fi(x*) 

0.009502 

0.009557 

0.020037 

f2(x*) 

0.001618 

0.0015899 

0.004813 

f3(x*) 

67.28 

66.801 

23.598 

0.009502 

0.0015899 

23.598 

i=l 
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Table  4.2 


Multiobjective  Optimization  of  Two-Bar  Truss 


Permissible  design  variable  change  =  -5% 

^  =  0.01,  C  =  10^.  x(o)  =  {ol 


Utility  Function  Method 

Lexicographic  Mediod 

Goal  Programming  Method 

Const.  Coef. 

Variable  Coef. 

Optimization  Order 

p=l 

P=2 

Optimal 

Design 

X,  =  0.1309 

0.1463 

0.15389 

0.056296 

0.14466 

0.051294 

0.051807 

Variables 

X* 

X2  =  0.12709 

0.1299 

0.13483 

0.056303 

0.1440 

0.051309 

0.051807 

f.(x*) 

0.01QS06 

0.009950 

0.009603 

0.018932 

0.009604 

0.020037 

0.019919 

f2(X*) 

0.001867 

0.001792 

0.001669 

0.004399 

0.001664 

0.004809 

0.004770 

f3(X') 

59.337 

63i49 

66.404 

25.898 

66.391 

23.599 

23.831 

iFi(X-) 

0.999438 

0.977291 

1.014075 

1.819317 

1.010398 

1.998782 

1.980792 
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Table  4.3 


Single  Objective  Optimization  of  Two-Bay  Truss 

Permissible  design  variable  change  =  -5% 

^  =  0.01,  Xi(o)  =  0.1  i=ltol0 


Minimization  of 

Objective 

PsT 

Ppr 

Weight.  W 

0.13816 

0.13765 

0.11274 

0.09648 

0.09339 

0.00100 

Optimal 

0.13782 

0.13777 

0.11315 

Design 

0.27661 

0.27637 

0.33788 

Variables 

0.10103 

0.09899 

0.00100 

x* 

0.27780 

0.27725 

0.33772 

i=l  to  10 

0.14537 

0.14671 

0.11810 

0.14417 

0.14507 

0.11884 

0.14808 

0.14998 

0.12110 

0.15119 

0.14920 

0.12122 

fi(X*) 

0.048694 

0.048701 

0.049354 

f2(X*) 

0.010019 

0.00981975 

0.0084567 

f3(X*) 

0.252557 

0.252307 

0.22730 
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Table  4.4 


Multiobjective  Optimization  of  Two-Bay  Truss 

Permissible  design  variable  change  =  -5% 

^  =  0.01,  Xi(o)  =  0.1  i=l,10 


Approach 

Utility  Function 
Method 

Lexicographic 

Method 

Goal 

Programming 

Method 

Const. 

Variable 

Optimization  Order 

p=2 

Coef. 

Coef. 

f1.f2.f3 

0.12247 

0.13753 

0.13697 

0.14197 

0.11389 

0.01779 

0.09207 

0.09588 

0.03602 

0.00100 

Optimal 

0.11406 

0.13603 

0.15071 

0.16162 

0.11142 

Design 

0.32120 

0.27907 

0.29918 

0.36978 

0.00100 

Variables 

0.01000 

0.09880 

0.07856 

0.00100 

0.00100 

xr 

0.33931 

0.27920 

0.29422 

0.34327 

0.33617 

i=l,10 

0.10835 

0.14554 

0.15018 

0.11642 

0.11907 

0.12838 

0.14278 

0.15061 

0.12668 

0.11901 

0.13045 

0.14917 

0.10706 

0.10488 

0.12277 

0.11834 

0.14972 

0.12839 

0.11305 

0.12420 

fl(X*) 

0.048769 

0.048742 

0.048788 

0.048664 

0.049360 

f2(X*) 

0.012034 

0.00997 

0.010307 

0.011721 

0.008479 

f3(X*) 

0.22903 

0.25196 

0.24841 

0.24438 

0.22736 

ZFiCX*) 

i=l 

2.431212 

1.946322 

2.094988 

2.694381 

0.999241 

m  =  1,0 

E  ■  10*10® 

Density  »  4.6 
Nominal  areas  =  0,1 


Figure  4.1  Two-Bar  Trass 
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Figure  4.2  SUbilUy  Robustness  Index  vs  the  Value  of  Design  Variables  for 
Two-Bar  Truss 
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Figure  4.10  Stability  and  Performance  Robustness  Indices  vs  the  Output 
Weighting  Matrix  Coefficient  for  Two>Bar  Truss 
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CHAPTER  5 


ACTUATOR/SENSOR  LOCATION  SELECTION 


The  problem  of  selection  of  optimal  locations  of  actuators/sensors  in  large  flexible 
structures  has  been  under  study  for  more  than  a  decade.  This  problem  was  considered  as  an 
integer  programming  problem  and  solved  using  heuristic  techniques  by  Haftka  and  Adelman 
(1985),  and  Skelton  and  DeLorenzo  (1983).  The  complexity  of  Skelton  and  DeLorenzo’s 
algorithm  for  m  actuators  on  n  available  positions  is  given  by 

Cos  =  Y(n+m)(n-m+l)  (5.1) 

The  complexities  of  WOBI  (worst-out-best-in)  and  ESPS  (exhaustive-single-  point- 
substitution)  algorithms,  studied  by  Haftka  and  Adelman,  can  be  expressed  as 

CwoBi  =  (5-2) 

and 

Cesps  =™^x(n-m)  (5.3) 

All  these  three  heuristic  techniques  yield  only  local  optimal  solutions  with  no  guarantee  about 
the  global  optimum. 

Recently,  Salama  et  al.  (1987)  considered  a  similar  problem  in  the  framework  of 
combinatorial  optimization  and  solved  it  using  the  simulated  annealing  heuristic  algorithm. 
The  conqjlexity  of  this  formulation  for  m  actuators  on  n  available  positions  can  be  expressed  as 

CcoM=C(n.m)=-j^j,  (5.4) 

Generally  speaking,  the  simulated  annealing  algorithm  will  converge  to  the  global  optimum, 
which  corresponds  to  the  lowest  energy  state,  if  an  appropriate  cooling  temperature  6),  is 
chosen.  The  selection  of  the  temperature,  B|,>  is  critical  in  the  simulated  annealing 
algorithm  for  which  there  is  no  systematic  approach.  Hence  a  trial  and  error  procedure  is 
usually  used  to  find  the  value  of  Oh- 

In  this  work,  three  additional  methods  are  investigated  for  solving  the  actuator  and  sensor 
location  problem. 

5.1  Problem  Formulation 

The  optimal  actuator/sensor  location  selection  problem,  by  nature,  is  a  zero-one  type 
combinatorial  optimization  problem.  The  system  equations  (2.1)  and  (2.2)  will  be  changed  if 
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the  locations  of  actuators/sensors  aie  changed.  Hence,  a  more  flexible  modeling  technique  is 
necessary  during  the  optimization  process.  For  consistency,  the  optimal  linear  quadratic 
regulator  is  used  for  the  control  system  design.  The  system  equations  and  the  control  gain  are 
expressed  in  Eqs  (2.2)  to  (2.6).  The  input  matrix  B  in  system  equation  (2.2)  can  be  expressed  as 


0 

M''D 


(5.5) 


in  physical  coordinates,  and  as 


(5.6) 


in  modal  space. 

The  inverse  of  the  input  weighting  matrix  R  is  assumed  to  be  a  diagonal  matrix  containing 
ones  and  zeros  only  with  I’s  corresponding  to  locations  with  actuators  and  O’s  corresponding  to 
those  without  actuators.  Hence,  the  system  equations  will  remain  unchanged  during  the 
optimization  process.  Similarly,  the  output  weighting  matrix  Q  can  be  modified  such  that  the 
system  equations  for  the  estimation  part  are  unchanged  during  the  optimization  process  for 
sensor  location  selection.  For  simplicity,  the  estimation  part  is  neglected,  and  only  the  actuator 
location  selection  problem  is  considered  in  this  work. 

The  objective  function  (criterion)  proposed  to  be  used  in  the  actuator/sensor  location 
selection  problem  is  the  energy  dissipated  by  the  active  controller  which  can  be  written  as 

Ec  =  |jq^Dcqdt  (5.7) 

^  0 


where  q  is  the  velocity  vector  and  D,.  is  the  induced  damping  matrix  by  the  active  controller. 
By  letting 


Pit  Pi2 
P21  P22 


(5.8) 


(5.9) 
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0  0 
~  [o  Dc 

or 

_  fO  0 

Dc  =  -  ,  -T 

0  P22 


(5.11) 


(5.12) 


(5.13) 


(5.14) 

and  P  is  the  solution  of  the  Lyapunov  equation 

Aci'*'P  +  PAci  =  Dc  (5.15) 

Equation  (5.12)  shows  that  the  energy  dissipation  depends  on  the  initial  state  which  may  not  be 
available  for  practical  problems.  A  simple  way  to  eliminate  the  dependence  on  the  initial  state 
is  to  average  it  out  by  assuming  the  initial  state  (qo  qo)  to  be  a  random  variable  uniformly 
distributed  on  the  surface  of  the  2n-dimensional  unit  sphere.  Hence,  an  upper  bound  on  the 
dissipation  energy  can  be  expressed  as 


the  dissipation  energy  can  be  rewritten  as 


^  /t 


T  T 

q  q 

J 

Wt 


Applying  the  transition  matrix  of  states  in  Eq.  (5.1 1)  gives 


E  =  — 


qo’^qo  j-le^^DcC^^dt-j 


qo 

qo 


which  can  be  simplified  as 


Ec  =  | 


T  •  ^ 

qo  qo 


qo 

qo 


where 


Aci  =  A  +  BG  =  A  -  BR-*B’*'P 


(5.10a) 


(5.10b) 
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(5.16) 


Let 


(5.17) 


be  the  inverse  of  the  input  weighting  matrix,  with  Fj  denoting  a  binary  variable  indicating  the 
presence  or  absence  of  an  actuator  at  position  i.  The  zero-one  optimization  problem  for  the 
actuator/sensor  location  selection  problem  can  then  be  expressed  as  follows: 

Maximize  Ec  (S.18) 

•r,ri  -  r.  • 

subject  to 


ri  +  r2  +  '  •  •  +  r„  =  m 


ri€(0,l),  i  =  1.2,...,n 


S2  Sequential'Best'Adding  (SBA)  Algorithm 

Sequential-Best-Adding  algorithm  is  a  simple  heuristic  search  procedure.  As  Skelton  and 
DeLorenzo’s  WOBI  and  ESPS  approaches,  the  sequential-best-adding  method  finds  a  local 
optimum.  However,  the  complexity  is  reduced  very  much,  especially  when  a  small  number  of 
actuators  (m)  are  used  in  a  large  feasible  space  (n). 

The  SBA  algorithm  can  be  described  as  follows: 

1.  Set  i=l,  j=0,  and  m=  specified  number  of  actuators  to  be  used. 

2.  Evaluate  the  objective  function  by  enumerating  all  possible  combinations  with  i 
actuators  where  j  (j^)  of  them  have  been  fixed  from  previous  calculations  in  the  feasible 
space. 

3.  Hnd  the  best  combination,  set  i=i+ 1  and  j=3+l . 

4.  If  i>m  then  stop  and  list  the  best  combination  as  the  solution,  otherwise  go  to  step  (2). 

This  approach  can  be  depicted  in  the  form  of  a  tree  structure  as  shown  in  Hg.  5.1  in  which 
only  one  branch  has  been  explored.  Hence,  the  complexity  of  the  SBA  algorithm  can  be 
expressed  as 
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CsBA  =  (5*  19) 

For  a  large  system,  the  parallel  processing  techniques  can  be  used  to  improve  the  efficiency  of 
the  sequential-best-adding  algorithm. 

5J  Penalty  Function  Method 

One  possible  way  to  solve  nonlinear  integer  programming  problems  is  use  a  penalty 
function  approach  in  which  a  penalty  is  included  for  non-integer  solutions.  This  relaxation 
technique  can  also  be  used  for  the  actuator/sensor  location  selection  problem.  The  relaxed 
version  of  equation  (5.18)  gives 

Maximize  Ec  +  Fp  (5.20) 

^.72  •••?•.} 


subject  to 


rj  +r2  +  •  •  •  +Fn  =  m 


where  Fp  is  the  penalty  function  because  of  the  non-zero-one  design  parameters. 

One  type  of  penalty  function  for  general  integer  or  discrete  optimization  problem  (Rao, 
1984)  is 


Fp=Yi 

i=sl 


xi-yi 


Zi-yi 


1- 


xj-yi 

Zi-yi 


P 


(5.21) 


where  xj,  (i=l,..,n)  are  design  parameters,  y  and  ^  (>1)  arc  constants  and  yi  and  zi  are  the  two 
neighboring  integer  values  of  Xi; 

yi  ^  Xi  ^  Zi  (5.22) 


For  a  zero-one  programming  problem,  yi  and  Zj  are  0  and  1,  respectively.  Hence  the  penalty 
function  in  Eq.  (5.21)  can  be  simplified  as 


4xi(l-Xi) 


(5.23) 
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This  penalty  function  is  shown  in  Fig.  5.2  for  different  values  of  y  and  p.  The  maximum 
magnitude  of  the  penalty  function  is  changed  when  the  value  of  y  is  changed,  while  its  shape  is 
adjusted  by  the  value  of  p.  The  choice  of  the  value  of  y  strongly  influences  the  convergence  of 
the  objective  function  while  that  of  p  does  not.  A  large  value  of  y  might  cause  to  find  only  a 
local  optimum.  Hence,  the  value  of  y  is  selected  small  enough  at  the  beginning  and  increased 
gradually.  The  value  of  p  is  a  constant  and  greater  than  1  in  order  to  maintain  the  continuity  of 
the  first  derivative  of  the  penalty  function  over  boundary  points. 


Another  form  of  penalty  function  is  an  exponential  form: 


Fp=YS 


i=l 


-a(l-x,) 


(5.24) 


This  function  gives  a  flatter  shape  penalty  function  on  non-feasible  solutions  as  shown  in 
Fig.  5.3.  As  with  the  penalty  function  of  Eq.(5.23),  the  value  of  y  in  Eq.(5.24)  strongly  affects 
the  convergence  of  the  optimization  process.  The  values  of  a  and  P  influence  the  shape  of  the 
penalty  function.  The  value  for  a  is  suggested  to  be  no  greater  than  10. 

5.4  Genetic  Algorithms 
5.4.1  Background 

Genetic  algorithms  (Davis,  1987  and  Goldberg,  1989)  are  basicly  guided  random  search 
techniques  derived  from  the  natural  genetics  of  populations.  The  decision  variables  (or  design 
parameters)  are  coded  as  a  string  of  binary  bits  which  correspond  to  the  chromosomes  in 
natural  genetics.  The  objective  function  value  corresponding  to  the  design  vector  plays  the  role 
of  fimess  in  natural  genetics.  The  artificial  recombination  among  the  population  of  strings  is 
based  on  the  fimess  and  the  accumulated  knowledge.  In  every  new  generation,  a  new  set  of 
strings  is  created  by  using  randomized  parents  selection  and  crossovers  from  the  old  set  of 
strings  (or  old  generation).  Although  randomized,  genetic  algorithms  are  not  siiiq)le  random 
search  techniques.  They  efficiently  explore  the  new  combinations  with  the  available  knowledge 
to  find  a  new  generation  with  better  fimess  or  objective  function  value. 

Genetic  algorithms  have  been  developed  and  studied  by  Holland  (1973)  and  his  students. 
These  algorithms  are  computationally  simple  but  powerful  in  their  search  for  improvement.  A 
growing  number  of  applications  can  be  found  in  science,  engineering,  business  and  social 
sciences. 

In  the  area  of  optimization,  three  basic  approaches  are  conunonly  used  for  solving 
nonlinear  programming  problems.  The  first  category  is  calculus  based  methods  which  use  the 
local  gradient  information  to  decide  the  search  direction  and  to  determine  the  convergence  of 
optimum.  These  approaches  find  local  optima  instead  of  global  ones  for  multimodal  functions. 
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They  are  not  useful  for  discrete  or  integer  programming  problems.  The  second  category  is 
enumerative  schemes  which  are  commonly  used  in  discrete  or  combinatorial  optimization. 
Although  these  methods  are  simple,  they  can  only  be  applied  to  a  problem  with  a  small  search 
space.  For  practical  problems,  the  search  space  usually  is  too  large  to  enumerate  the  entire 
space;  in  stead,  heuristic  algorithms  and  knowledge  based  search  are  used  to  increase  the 
search  efficiency.  The  third  category  belongs  to  random  search  techniques  which  have  attracted 
much  attention  when  the  previous  two  types  of  approaches  fail.  Random  search  can  be  divided 
into  unguided  and  guided  random  searches.  Unguided  random  searchs  can  be  expected  to  do  no 
better  than  enumerative  methods.  The  genetic  algorithm,  a  guided  random  search,  uses  random 
selection  as  a  tool  to  guide  a  highly  efficient  search  through  a  binary  coded  feasible  design 
space. 

The  genetic  algorithm  has  been  applied  to  different  optimization  problems  by  De  Jong 
(1980)  and  Goldberg  and  Samtani  (1986).  It  doesn’t  mean  that  genetic  algorithm  can  find  the 
optimal  solution  for  all  kinds  of  problems  efficiently.  On  the  other  hand,  genetic  algorithms 
usually  give  reasonable  near  optimal  solutions  (a  set  of  solutions).  Hence,  genetic  algorithms 
are  very  useful  to  solve  a  difficult  and  complex  problem  which  can  not  effectively  be  solved  by 
conventional  approaches.  Moreover,  since  the  whole  process  is  manipulated  on  a  set  of  points 
(population),  genetic  algorithms  are  suitable  for  parallel  processing  to  increase  the  efficiency. 

5.4.2  Genetic  Algorithms  in  Optimization 

A  simple  genetic  algorithm  involves  copying  strings  and  swapping  partial  strings  between 
two  mating  strings.  Three  basic  but  very  important  operators  -  reproduction,  crossover  and 
mutation  -  used  to  produce  new  generations  over  and  over  again. 

Reproduction  is  a  randonuzed  selection  process  in  which  individual  strings  are  copied 
according  to  their  fitness  (or  objective  function)  value.  Similar  to  the  natural  selection  process 
in  which  a  Darwinian  survival  of  the  fittest  among  string  creatures,  the  string  with  a  higher 
fitness  value  has  a  higher  probability  to  be  reproduced.  The  probability  of  reproduction  can  be 
determined  by  dividing  the  individual  fitness  by  the  sum  of  fimesses  of  the  ciurent  generation. 
For  example,  if  a  set  of  strings  with  a  population  size  5  have  fimesses  SO,  75,  100,  125,  and 
150,  their  probabilities  of  reproduction  are  50/5(X),  75/5(X),  100/500,  125/500,  and  150/500, 
respectively.  (Dnce  a  string  is  chosen  for  reproduction,  a  duplication  is  made  and  the  replica  is 
sent  to  the  mating  pool  for  further  genetic  operator  actions. 

Crossover  (or  recombination)  is  a  primary  operator  in  the  mating  process  which  generates 
the  offsprings  or  new  generations.  Two  steps  are  involved  in  the  crossover  process.  The  first 
step  is  randomly  selecting  the  position  to  break  (crossover)  between  two  mating  couple.  The 
second  step  is  reunion  of  these  two  mating  couple.  In  natural  genetic  theory,  the  recombination 
fiaction  depends  on  the  distance  between  the  chromosomes  of  the  mating  couple.  The  closer 
the  two  loci  are,  the  less  likely  that  such  a  breakage  and  reunion  will  occur  between  them;  the 
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farther  apart  they  are,  the  more  likely  such  a  recombination  will  occur.  For  the  simple  artificial 
genetic  algorithm,  a  fixed  and  preassigned  probability  of  crossover  is  used.  For  example,  let 
two  strings  A)  and  A2  with  8-bit  be  selected  and  reproduced  in  the  mating  pool  as: 

Aj:  01100100 
A2: 10011100 

Seven  crossover  positions  are  possible  to  select.  If  the  crossover  point  is  chosen  at  position  4 
by  a  random  choice,  the  two  new  strings  A'l  and  A'2  are  generated  to  be  a  part  of  new 
generation  as: 

Ai:  011010100  A'pOllOllOO 

A2: 100111100  A'2: 10010100 

where  |  in  Ai  and  A2  denotes  the  crossover  posidon.  This  process  can  be  implemented  by  (i) 
selecting  the  crossover  posidon  by  a  uniformly  distributed  random  process  between  1  and  the 
string  length  minus  1  ,  i.e.,  [1,1-1]  and  (ii)  swapping  the  contents  of  mating  couple  after  the 
crossover  point. 

Mutadon  plays  a  secondary  role  in  the  operadon  of  genedc  algorithms.  Mutadon,  changing 
a  particular  bit  of  coded  string  from  0  to  1  or  vice  versa,  is  a  random  walk  through  the  string 
space.  In  natural  genedcs,  mutadon  is  the  ultimate  source  of  genedc  variation  which  can  arise 
only  if  there  is  some  change  in  the  genedc  material.  Because  spontaneous  mutadon  rates  are 
typically  quite  small  (on  the  order  of  KT*  to  lOr^  mutations  per  locus  per  generation),  the 
effect  is  very  small  over  the  course  of  a  few  generations.  On  the  other  hand,  the  cumulative 
effects  of  mutation  over  a  long  period  of  time  can  become  appreciable.  Hence,  mutation  is 
needed  in  artificial  genetic  systems  to  act  as  an  insurance  policy  against  premature  loss  of 
important  information. 

The  genetic  algorithm  can  be  applied  for  solving  optimization  problems  as  follows: 

1.  An  appropriate  chromosome  representation  should  be  defined  to  represent  the 
combinations  of  design  parameters  which  correspond  to  the  fitness  or  objective  function 
values.  The  representation  should  be  one-to-one  mapping  in  order  to  have  a  normal 
coding  and  decoding  processes. 

2.  The  population  size  and  the  maximum  number  of  generations  should  be  specified.  The 
probabilities  of  crossover  and  mutation  are  selected.  A  set  of  initial  population  in  the 
genetic  system  will  also  be  generated. 

3.  Evaluate  the  fimess  (or  objective  function)  value  of  each  individual  in  the  current 
generation.  The  objective  function  plays  the  role  of  the  environment  to  decide  the  fimess 
of  a  chromosome.  Two  important  steps  are  usually  taken  to  modify  the  objective 
function  values  for  fimess  evaluation.  The  first  step  is  scaling  to  change  the  distribution 
range  of  the  fimess  value.  The  second  one  is  normalization  process  which  can  change  the 
speed  of  fixation  of  population.  For  this  reason,  users  of  genetic  algorithms  have 
employed  a  variety  of  normalization  techniques.  If  normalization  process  stresses 
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improvements  too  much  it  will  lead  to  a  premature  dominance  of  a  single  string.  When 
this  happens  crossover  becomes  of  little  value,  and  the  algorithm  ends  up  in  intensively 
searching  the  solution  space  in  the  region  of  the  last  good  individuals  found.  If  it  doesn’t 
stress  good  fimess,  the  algorithm  might  fail  to  converge  to  a  good  result  and  will  be  more 
likely  to  lose  the  best  combination  of  its  population.  Hence,  the  normalization  process  is 
critical  in  genetic  algorithms. 

4.  Apply  the  three  operators  •  reproduction,  crossover  and  mutation  -  on  the  old  generation 
to  generate  the  new  population  for  the  next  generation. 

5.  Repeat  step  (3)  and  (4)  until  the  maximum  number  of  generations  is  achieved. 

5.4.3  Some  Probability  Relations  in  Genetic  Algorithms 

To  understand  why  the  genetic  algorithms  work,  some  probability  study  would  be  helpful. 
Three  operators  without  restriction  are  considered  in  this  section.  A  schema  in  genetic 
algorithms,  chromosomes  with  specified  0-1  combination,  is  defined  as  the  symbol  set  {0,1,*) 
in  which  ♦  denotes  ’don’t  care’.  To  quantitatively  evaluate  the  probability  of  survival  of  a 
specific  schema,  the  schema  order  and  defining  length  are  needed.  The  order  of  a  schema  A, 
o(A),  is  the  number  of  fixed  positions  (i.e.  O’s  and  I’s)  present  in  the  template.  For  example, 
Al=011*0***,  A2=l*l****l  have  schema  orders  o(Al)=4  and  o(A2)=3,  respectively.  The 
defining  length  of  a  schema  A,  S(A).  is  the  distance  between  the  first  and  the  last  specific  string 
positions.  For  example,  A1  and  A2  have  the  defining  lengths  5(A1)=5-1=4,  and  5(A2)=8-1=7, 
respectively. 

During  reproduction,  a  string  is  selected  with  the  probability  of  its  own  fimess  divided  by 
the  total  fimess  in  that  generation.  After  picking  n  individuals  firom  the  population  which  have 
the  the  probability  P(A,i)  in  which  A  indicates  the  specific  schema  and  i  denotes  the  i*^ 
generation,  P(A,i+l)  can  be  expressed  as  follows: 

P(A.i+l)  =  P(A,i)xn-^  (5.25a) 

I 

or 

P(A,i+l)  =  P(A,i)-^  (5.25b) 

f 

where  f(A)  and  f  denote  the  fimess  of  specific  schema  and  average  fimess,  respectively. 

During  crossover,  a  schema  can  survive  if  the  crossover  site  falls  outside  the  defining 
length.  Hence  the  survival  probability  under  simple  crossover  is  given  by 
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(5.26) 


where  1  is  the  length  of  the  string  A.  If  a  random  choice  is  used  in  crossover  and  the 
probability  of  crossover  is  selected  as  Pc.  the  survival  probability  is  given  by  the  expression 


Ps^l-Pc 


S(A) 

1-1 


(5.27) 


Muution  is  the  random  alternation  of  a  single  position  with  a  specified  probability  P^.  A 
specific  schema  A  can  survive  if  and  only  if  all  alleles  (elements)  at  the  specified  positions 
must  themselves  survive.  A  single  allele  survives  with  a  probability  of  1-Pn,  and  hence  the 
survival  probability  of  schema  A  because  of  mutation  is 

(5.28) 


Since  Pm  is  usually  very  small  compared  to  1,  Eq.  (5.28)  can  be  simplified  as 

P«n  =  (l-o(A)xP„) 


(5.29) 


Therefore,  the  probability  of  a  specific  schema  in  the  next  generation,  under  three  genetic 
operators,  can  be  obtained  by  combining  equations  (5.25),  (5.27),  and  (5.29): 

f 


P(A,i+l)  ^  P(A,i)^^ 
f 


1-R 


8(A) 

1-1 


[l^(A)xP„j 


(5.30) 


or 


P(A,i+l)^P(A,i)-^i^ 


(5.31) 


Since  P^  is  very  small  compared  to  1  or  Pc,  and  1,  f(A),  Pc,  and  P^  are  constants,  6(A) 
becomes  much  more  important  than  o(A).  Hence,  the  fitness  and  the  defining  length  of  string  A 
become  the  two  factors  to  decide  the  survival  probability.  An  observation  of  Eq.(5.31)  shows 
that  the  short,  low-order,  and  above-average  schemata  receive  an  exponentially  increasing 
probability  in  subsequent  generations. 

5.4.4  More  on  Simple  Genetic  Algorithms 

For  the  actuator/sensor  location  selection  problem,  a  type  of  combinatorial  optimization 
problem,  the  genetic  algorithm  will  be  a  suitable  approach.  The  chromosomal  representation  of 
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the  design  vector  simply  involves  using  a  string  of  binary  bits.  The  locations  in  the  binary 
string  indicate  the  numbered  positions  on  the  controlled  structure  to  place  the  actuators  or 
sensors.  The  I’s  and  O’s  in  the  binary  string  denote  the  presence  and  absence  of  the 
actuators/sensors.  For  example,  the  string  10100100  means  that  three  actuators  (or  sensors)  are 
chosen  and  located  at  positions  1,3,  and  6  out  of  the  possible  8  positions. 


Several  modifications  are  made  to  the  genetic  algorithm  described  in  section  5.4.2  before 
using  it  for  the  solution  of  the  actuator/sensor  location  selection  problem: 

1.  The  initial  population  is  randomly  generated  with  the  restriction  that  no  two  individuals 
be  allowed  to  have  the  same  chromosome  in  order  to  get  the  maximal  variety. 

2.  The  successive  flips  of  a  biased  coin  were  used  to  generate  the  initial  population 
efficiently. 

3.  A  normalization  and  scaling  process  has  been  applied  to  increase  the  differences  between 
the  fitnesses  of  different  chromosomes.  The  normalization  process  is  indicated  by  the 
relation: 


f'(x)  = 


^mix  f(x) 
f  -f  - 

^mtx  ^mm 


(5.32) 


where  f(x)  is  the  objective  function  value  with  design  vector  x  and  f  (x)  is  the  normalized 
fitness  value  for  chromosome  x. 

4.  A  constraint  is  placed  on  the  total  number  of  actuators  used.  Two  approaches  can  be  used 
for  this  purpose.  In  the  first  approach,  a  restriction  is  placed  on  the  crossover  process 
which  is  permitted  to  generate  two  children  with  the  specified  number  of  actuators  (i.e., 
number  of  I’s  in  the  chromosome).  Otherwise,  the  crossover  process  is  not  allowed.  In 
the  second  approach,  no  restriction  is  placed  on  the  crossover  process;  but  a  corrective 
action  is  taken  afterwards.  After  mating,  the  number  of  I’s  in  the  two  children  are  made 
equal  to  the  specified  number  by  a  random  mutation  process.  For  example,  let  two 
mating  couple  with  chromosomes  be  given  by 

Aj:  10100010 
A2:00101010 

where  the  total  numbo’  of  I’s  is  specified  as  3.  According  to  the  first  method,  the  only 
possible  positions  for  crossover  are  5  (i.e.,  the  position  between  5*  and  6*  alleles  from 
left  side),  6,  and  7;  and  none  of  them  gives  a  new  chromosome  (or  combination).  On  the 
other  hand,  no  more  improvement  can  be  made  through  these  two  mating  couple. 

If  the  second  approach  is  chosen,  crossover  can  be  made  at  any  position  among  1  to  7 
with  a  probability  y .  If ,  for  example,  position  4  is  randomly  chosen  for  crossover,  then 
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Ai:  10100010  P=y  A'l:  10101010 


(5.33) 


A2: 00101010  A'2: 00100010 


The  numbers  of  I’s  in  A'l  and  A'2  do  not  satisfy  the  constraint  (value  of  3)  and  hence  a 
random  mutation  should  be  used  on  A'l  (A'2)  to  reduce  a  1  in  A'l  (increase  a  1  in  A'2). 
Suppose  the  optimal  combination  is  "10101000",  then  the  probability  of  generating  the 


optimal  offspring  is 


28 


A'l:  10101010 

4 

10101000 

(5.34) 

A'2: 00100010 

^^4 

10100010 

Ai:  10100010 

A2: 00101010 


(5.35) 


It  is  obvious  that  the  second  method  gives  more  freedom  in  the  crossover  process  and 
will  reduce  the  probability  of  non-improving  crossover. 

5.  A  variable  population  size  is  used  in  the  actuator/sensor  location  selection  problem. 
Since  the  genetic  algorithm  is  highly  dependent  on  the  normalization  technique  used,  the 
best  chromosome  which  was  found  since  the  beginning  of  the  process  is  not  guaranteed 
to  be  included  in  the  cuirent  generation.  To  preserve  the  best  individual,  the  best 
individual  in  the  old  generation  (A)  is  duplicated  to  new  generation  if  there  is  no  other 
individual  who  is  at  least  as  good  as  A.  Hence,  the  population  size  is  increased  by  1  if  a 
duplication  is  made. 

6.  Although  random  mating  is  the  most  important  mating  system  in  many  natural 
populations,  there  are  certain  departures  from  the  random  mating  that  can  also  be  very 
important.  The  definitions  and  characteristics  of  several  important  mating  systems  are 
listed  below; 

Random  mating:  Choice  of  mates  independent  of  genotype  and  phenotype. 

Positive  assortative  mating:  Mates  phenotypically  more  similar  than  would  be  expected 
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by  chance. 

Negative  assortative  mating:  Mates  phenotypically  more  dissimilar  than  would  be 
expected  by  chance. 

Inbreeding:  Mating  between  relatives. 

In  humans,  positive  assortative  mating  occurs  for  height,  IQ  score  and  certain  other 
traits.  Negative  assortative  mating  is  apparently  quite  rare,  but  an  obvious  example  is 
sex.  The  principal  effect  of  inbreeding  in  a  population  is  to  increase  the  frequency  of 
homozygous  genotypes  at  the  expense  of  the  frequency  of  heterozygous  genotypes.  This 
effect  can  most  easily  be  seen  in  the  case  of  repeated  self-fertilization.  Detrimental 
effects  of  inbreeding  (called  inbreeding  depression)  are  found  in  virtually  all  outcrossing  * 

species,  and  the  more  intense  the  breeding,  the  more  harmful  the  effects.  The  inbreeding 
might  happen  in  random  mating  process  with  no  restriction.  The  frequency  (or 
probability)  of  inbreeding  will  increase  when  the  population  with  the  same  genotype 
increases.  In  artificial  genetic  algorithms,  the  inbreeding  process  gives  no  new 
information  and  reduces  the  time  to  fixation.  Hence,  the  restriction  on  inbreeding  should 
be  considered. 

5.5  Examples 

To  illustrate  the  feasibility  of  the  three  methods  discussed  in  the  previous  sections,  a 
two-bay  truss  (Fig.  2.3)  and  a  six-bay  truss  (Fig.  2.9)  are  considered  as  exanq)les.  The 
number  of  actuators  is  assumed  to  be  3.  The  optimal  linear  quadratic  regulator  is  applied 
to  solve  the  optimal  control  gain.  The  dissipation  energy  of  active  controller  is  used  as 
the  objective  function  for  maximization.  The  output  weighting  matrix  is  assumed  to  be 
1(X)0*I  (I  is  the  identity  matrix)  in  both  the  examples.  The  input  weighting  matrix  has 
fixed  dimension  of  n  (n  is  the  total  number  of  d.o.f.’s  in  the  structure)  and  is  composed  as 
described  in  section  5.1. 

By  using  the  sequential-best-adding  algorithm,  the  result  of  the  two-bay  truss  is 
presented  in  Table  5.1.  The  optimal  solution  for  level  1  (i.e.,  when  one  actuator  is  used) 
is  at  position  2.  The  optimal  solution  far  level  2  (when  two  actuators  are  used)  is  at 
positions  2  and  4  which  is  a  child  in  branch  2.  But  the  optimal  solution  for  level  3  (when 
three  actuators  are  used)  is  at  positions  2,  4,  and  7  which  is  not  same  as  the  global 
optimal  solution  (1,3,5).  The  number  of  function  evaluations  used  is  21  out  of  a  total 
combination  of  56  (  c(8,3)^!/3!5! ).  For  the  six-bay  truss,  the  global  optimal  solution 
for  3  actuators  is  at  (2,3,4).  The  SBA  approach  gives  the  optimal  solutions  for  levels  1  to 
3  at  (2),  (2,4),  and  (2,3,4),  respectively.  The  number  of  function  evaluations  used  is  69  ( 
=24+23+22)  out  of  a  total  combination  of  2024  ( c(24,3)=24!/21!3!). 

The  optimal  selections  of  actuator  locations  for  the  two-bay  truss  using  the  penalty 
function  method  with  different  parameters  are  shown  in  Table  5.2.  None  of  the 
parameters  give  the  global  optimal  solution  (1,3,5).  Also,  the  number  of  function 
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evaluations  is  much  higher  than  the  total  number  of  combinations  (S6).  Hence,  this  is  not 
a  feasible  approach  for  the  actuator/sensor  location  selection  problem. 

The  genetic  algorithm  is  also  applied  to  the  same  examples.  For  the  two-bay  miss, 
the  initial  population  size  and  the  maximum  number  of  generations  are  selected  as  8  and 
15,  respectively.  Results  obtained  with  different  crossover  probability  (from  0.6  to  1.0) 
and  mutation  probability  (from  0.001  to  0.005)  are  shown  in  Table  5.3.  The  number  of 
function  evaluations  in  10  generations  varies  from  18  to  35.  Most  of  them  find  the  global 
optimal  solution  (1,3,5).  A  small  number  of  them  find  the  second  optimal  solution 
(1,3,7)  which  is  in  the  same  branch  (1,3)  of  the  global  optimal  solution.  Three  randomly 
selected  initial  generations  are  used  for  comparison  in  two  cases;  (i)  Pc  =0.7  and 
Pm=0.001,  (ii)  Pc=0.8  and  Pm=0.001.  The  results  are  shown  in  Figs.  5.4  to  5.6  for  case 

(i)  and  in  Figs.  5.7  to  5.9  for  case  (ii).  Figure  5.4  shows  that  the  total  number  of 
genotypes  (configurations)  generated  is  constant  after  generation  32,  i.e.,  no  new 
genotypes  are  generated.  Hgures  5.5  and  5.6  show  that  the  average  and  the  maximal 
fitnesses  increase  with  the  number  of  generations  until  the  fixation  occurs.  The  same 
trend  is  revealed  in  Figs.  5.7  to  5.9  also. 

For  the  six-bay  truss,  the  initial  population  size  and  the  maximum  number  of 
generations  are  chosen  as  20  and  40,  respectively.  The  results  obtained  with  different 
crossover  and  mutation  probabilities  are  shown  in  Table  5.4.  The  number  of  function 
evaluations  in  15  generations  varies  from  110  to  271.  Most  of  them  find  the  global 
optimal  solution  with  the  acmators  at  positions  (2,3,4).  Few  of  them  find  the  second 
optimal  solution  at  (1,2,4)  which  is  in  the  same  branch  (2,4)  of  the  global  optimal 
solution.  Three  randomly  selected  initial  generations  are  used  for  comparison  in  two 
cases:  (i)  Pc=0.7  and  Pnj=0.001,  and  (ii)  Pc=0.8  and  Pn,  =0.001.  The  relations  between  the 
total  number  of  genotypes,  average  fimess,  maximal  fimess  and  the  number  of 
generations  are  shown  in  Figs.  5,10  to  5;12  for  case  (i)  and  in  Figs.  5.13  to  5.15  for  case 

(ii) .  The  premature  phenomenon  occurs  in  case  (i)  (line  c  in  Hgs.  5.10  to  5.12).  The  only 
way  to  reduce  the  possibility  of  prematurity  is  to  reduce  the  probability  of  inbreeding. 

5.6  Concluding  Remarks 

Three  approaches  are  presented  for  solving  the  actuator/sensor  location  selection 
problem.  The  sequential-best-adding  method,  although  no  global  optimal  solution  is 
guaranteed,  reduces  the  computational  effon  to  a  large  extent.  The  complexity  is 
proportional  only  to  the  number  of  variables  involved  in  the  optimization  problem.  For 
the  two-bay  truss,  SBA  gives  local  optimal  solution  which  is  the  third  global  optimal 
solution,  while  SBA  gives  the  global  optimum  for  the  six-bay  truss.  For  a  large  flexible 
structure  with  a  small  number  of  actuators  or  sensors,  the  SBA  approach  will  be  very 
efficient. 
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The  penalty  function  method  complicates  the  zero-one  actuator/sensor  location 
selection  problem  and  cannot  converge  to  the  best  solution  or  local  optimum. 

The  genetic  algorithm,  although  a  simple  guided  random  search  algorithm,  is  very 
effective  in  solving  actuator/sensor  location  selection  problem  when  posed  as  a  zero-one 
programming  problem.  In  most  cases  of  the  examples  considered,  the  global  optimal 
solution  is  found  in  reasonable  number  of  function  evaluations.  The  following 
recommendations  can  be  made  for  the  genetic  algorithm  (GA): 

1.  The  size  of  population  should  be  close  to  the  length  of  the  chromosome  used  in  the 
genetic  algorithm. 

2.  The  number  of  generations  needed  in  GA,  in  general,  is  a  state  of  art.  An 
interactive  procedure  to  decide  the  termination  of  GA  will  be  more  profitable. 

3.  The  probability  of  crossover  should  be  high  and  the  that  of  mutation  should  be 
low.  The  suggested  values  for  Pc  and  are  0.7  to  0.8  and  0.001  to  0.005, 
respectively. 

4.  No  specific  convergence  criterion  has  been  used  in  GA.  When  the  best  individual 
with  the  maximal  fimess  remains  unchanged  for  several  generations,  say  5  or  10, 
we  can  assume  that  the  process  has  converged. 

5.  As  shown  in  the  examples,  GA  gives  very  good  results.  However,  in  general,  GAs 
do  not  guarantee  convergence  to  the  global  optimal  solution.  A  hybrid  scheme  is 
suggested  for  solving  the  actuator/sensor  location  selection  problem.  The  hybrid 
procedure  to  be  followed  as:  (i)  start  with  the  GA  and  find  the  best  individuals  in 
the  last  generation,  (ii)  explore  the  branches  found  in  (i)  as  in  SBA  approach,  and 
(iii)  choose  the  best  local  optimum  fix>m  step  (ii). 

6.  To  improve  the  efficiency  of  the  procedure,  parallel  processing  techniques  can  be 
used  in  genetic  algorithms. 
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Table  5.1  Objective  Function  Velues  for  Two-Bay  Tnus 


Level  1  (•10') 

Levd2(«10^) 

Level  3  (*10*) 

(U) 

1.9S3S 

(144) 

1.9493 

(14,4) 

1.896S 

(144) 

1.8344 

(144) 

1.7293 

(14,7) 

1.7630 

(14,8) 

1.7138 

(1,3) 

2.3262 

(14.4) 

1.9409 

(144) 

11149 

(1.3.6) 

1.9114 

GESl 

(1.3.8) 

1.9207 

(1.4) 

1.9370 

(1.44) 

1.908S 

(1.4.6) 

1.6S6S 

(1.4.7) 

1.7570 

(W.8) 

1.6312 

(D* 

2.1469** 

(14) 

2.1020 

(144) 

I.6S99 

(14.7) 

1.8428 

(144) 

1.6988 

(1.6) 

1.8417 

(1.6.7) 

1.6166 

(1.6.8) 

1.8432 

(1.7) 

1.9694 

(1.7.8) 

1.6402 

(1.8) 

1.8973 

(24) 

1.9310 

(24.4) 

1.8899 

(244) 

1.7509 

(24.6) 

1.6888 

(24.7) 

1.8433 

(244) 

1.6663 

(2.4) 

3.2917 

a44) 

1.9393 

04.6) 

1.7231 

04.7) 

1.9642 

(2.4,8) 

1.7123 

(2) 

6.78S5 

(24) 

1.8886 

(24.6) 

14272 

03.7) 

1.6079 

(244) 

14622 

(24) 

1.8492 

06,7) 

14111 

(244) 

1.7079 

(2.7) 

1.9799 

(2.7.8) 

14303 

(2.8) 

1.8328 

(3.4) 

1.9S03 

(3.44) 

14044 

WESS 

(3,4,7) 

1.8440 

(3.4.8) 

1.6748 

(3) 

2.1173 

(34) 

2.0610 

044) 

1.6169 

(34,7) 

1.8993 

(344) 

1.6343 

Note:  *  -  Actuator  locations 

**  -  The  value  of  objective  function 
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Table  5.1  Objective  Function  Values  for  Two-Bay  Truss  (Continued) 


Level2(*10^)  j 

Uvel  3  (*10^) 

(3) 

2.1173 

(3.6) 

1.8599  1 

(3.6.7) 

1.6556 

(3.6.8) 

1.8320 

Bi 

lEuElHi! 

(3.7.8) 

1.6832 

B 

WMSMM 

Bi 

(4) 

4.0737 

(4.5) 

1.9728 

(4.5.6) 

1J046 

(4.5.7) 

1.6698 

(4.5.8) 

1.5336 

(4.6) 

1.8051 

(4.6,7) 

1J231 

(4,6.8) 

1.6609 

(4.7) 

1.9847 

(4,7.8) 

1J365 

(4.8) 

1.7925 

(5.6) 

1.6224 

(5.6,7) 

1.4662 

(5.6,8) 

1.6449 

(5) 

2.1772 

(5,7.8) 

1.5551 

(6) 

6,0217 

(6,7) 

1.6050 

(6,7.8) 

1.6324 

(6,8) 

2J064 

(7) 

2.2810 

(7.8) 

1.6750 

(8) 

6.1063 

1 
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Table  5.2  Penalty  Function  Method  for  Two^Bay  Truss 


Penalty  function  1 

Penalty  fiuction  2 

y 

0.2 

0.02 

0.0512 

0.131 

0.2 

0.819 

0.32 

2.2 

2.2 

1.408 

2.2 

2.2 

2.2 

■sa 

0 

0 

0 

0 

1 

0 

lEH 

0 

0 

0 

■■ 

0 

1 

0 

0 

0 

0 

■a 

0 

0 

Xi 

1 

1 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

1 

1 

1 

0 

1 

mm 

1 

0 

1 

1 

0 

0 

■9 

1 

1 

1 

1 

1 

1 

1 

f(x) 

(•10*) 

1.54 

1.662 

1.633 

nm 

1.651 

1.661 

1.562 

number 
of  eva¬ 
luations 

443 

292 

359 

367 

369 

363 

364 
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Tible  5.3  Genetic  Algorithm  for  Two^Baj  Tnue 


H 

m 

Best  final 
genotype 

Total  #  of 
evaluations 

Ofaj.  function 

max.  fiuwss 
found 

1.0 

TiTTirrrrM 

28 

2.114964-5 

1.0 

0.002 

10101000 

22 

2.114964-5 

1.0 

0.005 

10101000 

31 

2.1149C4-5 

1.0 

0.9 

0.001 

10101000 

35 

2.1149C4-5 

0.002 

10101000 

35 

2.1149C4-5 

1.0 

0.005 

10101000 

35 

2.114964-5 

1.0 

0.8 

0.001 

27 

2.065464-5 

0.9335 

0.002 

10100010 

28 

2.065464-5 

0.9335 

10101000 

32 

2.114964-5 

1.0 

0.7 

0.001 

10101000 

25 

2.114964-5 

1.0 

0.002 

10101000 

26 

2.114964-5 

1.0 

10100010 

32 

2.065464-5 

0.9335 

0.6 

0.001 

10101000 

23 

2.114964-5 

1.0 

0.002 

18 

1.845264-5 

0.6659 

0.005 

10101000 

34 

2.114964-5 

1.0 

Note:  *  ~  global  optimum,  **  ~  second  global  optimum 
-  local  optimum 
Populadou  size  »  8 
Number  of  generations  s  is 
Total  number  of  genotypes  (combinations)  =  56 
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Table  6.4  Genetic  Algorithm  for  Six-Bay  Tnus 


m 

m 

Best  final 
genotype 

Total  #  of 
evaluations 

Obj.  function 

max.  fitness 
found 

1.0 

0.001 

189 

9.372e+4 

1.0 

0.002 

I 

269 

9.372C+4 

1.0 

0.005 

I 

271 

9.372C-I4 

1.0 

0.9 

0.001 

I 

203 

9.372C-H4 

1.0 

0.002 

genotype  Il“"' 

202 

9.204044 

0.9586 

I 

271 

9.372e44 

1.0 

0.8 

0.001 

I 

142 

9.3720+4 

1.0 

0.002 

I 

125 

9.3720+4 

1.0 

0.005 

I 

245 

9.372e+4 

1.0 

0.7 

0.001 

m^Qniiiiiiii 

110 

9.204O+4 

0.9586 

0.002 

I 

141 

9.372C+4 

1.0 

0.005 

nmnTjiiiin 

194 

9.204O+4 

0.9586 

0.6 

0.001 

I 

177 

9.3720+4 

1.0 

0.002 

I 

205 

9.3720+4 

1.0 

0.005 

I 

198 

9.3720+4 

1.0 

**  genotype  n  a  1011 
Population  size  »  20 
Number  of  generatims  •  40 
Total  number  of  genotypes  (combinations)  =  2024 


;  second  optimum 
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SBA  algorithm  for 

m-actuator  from  n-feasible  locations 


•  • 

Lavtl  m 


^ ••  Sn) 


best  in  level  m 


Figure  5.1  Sequentiei-Beet-Adding  Algorithm  Search  Tree  Structure 
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Generation 


Accumulated  Genotypes  (Eraluations)  for  2-Ba7  Truss  with 
Different  LC.’s  (Pe=0.7,  P-=0.001) 


1- 


0“^ -S':"**  T™.  wltiDMfc^t  LC 
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N 


Fipxre  6.7 


Accumulated  Genotypes  (Evaluations)  for  2-Bay  Truss  with 
Different  LC.’s  (Pe=0.8,  Pni=0.00l) 


141 


2 


Figure  5.10  Accumulated  GenotTpes  (Evaluations)  for  O-Bay  Truss  with 
Different  LC.'s  (P^sO.?,  Pg, =0.001) 
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Figvra  5.12  Maximal  Fitness  Change  for  6>Ba7  Truss  with  Different  I.C.’s 
{Pe=0.7,  Pn=0.001) 


46 


Figure  5.13  Accumulated  Genotypes  (Eraluations)  for  6-Bay  Truss  with 
Different  LC.’s  (Pc=0.8,  =0.001) 
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Flfure  5.15  Maximal  Fitness  Change  for  S-Bay  Trass  with  Different  I.C.’s 
(Pe=0.8,  Pn=0.001) 
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CHAPTER  6 


CONCLUSIONS 


6.1  Work  in  Retrospect 

The  stnicmral  and  control  design  problem  of  flexible  structures  is  studied.  The 
substructures  concept,  widely  used  in  large/complex  structural  analysis,  is  used  to  fonnulate 
the  control  system  model  so  that  the  order  of  system  model  for  each  substructure  is  small 
enough  for  a  controller  design  with  no  further  model  reduction  needed.  Consequently,  the 
interacting  substructure  decentralized  control  procedure  will  give  a  suboptimal  control  gain 
without  the  spillover  problem.  The  interacting  force  that  is  due  to  the  active  controller  is 
transmitted  through  the  boundary  nodes  to  the  connected  substructures  and,  hence,  a  modified 
system  model  for  connected  substructures  will  be  made  for  control  design.  An  iterative 
procedure  is  performed  for  finding  the  control  gain  matrix  for  each  substructure.  Appropriate 
assembly  of  control  gains  for  substructures  gives  the  global  control  gain.  This  procedure 
produces  a  higher  gain  compared  with  that  of  LQR  designed  for  a  complete  system  model.  A 
unified  passive  damper  design  procedure  leads  to  an  approximation  in  the  higher  order  modes 
without  a  significant  change  in  the  lower  order  modes.  When  the  dual  passive/active  control 
(UPD/ISDC)  is  used,  the  system  reliability  and  lobusmess  are  increased.  The  performance 
index  is  also  reduced  by  a  large  amount  depending  on  the  quantity  of  passive  damping  added. 

Since  uncertainties  are  involved  in  the  model  for  control  system  design,  the  robustness  of 
stability  and  performance  should  be  considered  for  the  closed-loop  system.  The  stability  and 
performance  robustness  indices  are  defined  and  multiobjective  optimization  techniques  are 
introduced  to  solve  the  integrated  structural/control  design  problem.  The  major  motivation  of 
combined  structural/control  design  is  that  the  best  control  design  based  on  an  existing  structure 
might  not  constitute  the  best  closed-loop  system.  This  complicates  the  design  procedure,  but 
will  reduce  the  number  of  design  iterations  and  leads  to  a  better  design.  The  utility  function, 
lexicographic  and  goal  programming  methods  are  applied  to  solve  the  multiobjective  nonlinear 
programming  problem.  No  unique  optimal  solution  is  obtained  because  of  the  multimodal 
nature  of  the  objective  function. 

The  location  selection  problem  for  sensors/actuators  in  large  flexible  structures  is  very 
important  and  contains  a  large  number  of  possible  combinations  for  a  large  system.  Several 
heuristic  approaches  have  been  developed  by  previous  researchers.  The  sequential-best-adding 
algorithm,  penalty  function  method,  and  genetic  algorithm  are  introduced  in  this  woik.  The 
sequential-best-adding  approach  reduces  the  computational  effort  to  a  minimum  but  yields  only 
a  local  optimum.  The  penalty  function  method  is  not  appropriate  for  this  problem.  The  genetic 
algorithm  seems  to  be  the  most  promising  approach  which  is  demonstrated  on  simple 
examples.  The  genetic  algorithm  is  very  simple  to  implement  and  a  near  global  optimal 


solution  is  always  assured. 

6.2  Recommendation  for  Future  Research 

(1)  Investigate  the  experimental  implementation  of  interacting  substructure  decentralized 
control  of  large  flexible  structures. 

(2)  For  integrated  structural/control  design  problem,  investigate  the  feasibility  of  using  the 
genetic  algorithm  in  multiobjective  optimization  can  be  investigated. 

(3)  More  studies  on  genetic  algorithms  are  needed  for  application  to  the  sensor/actuator 
location  selection  problem  of  very  large  structural  systems.  A  parallel  processing  approach  in 
genetic  algorithms  is  also  worthy  of  investigation. 
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